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INFLUENCE OF LARGE DEFORMATIONS AND MIDPLANE FORCES 


ON THE PLASTIC BEHAVIOR OF IMPULSIVELY 
LOADED PLATES 
By 

Edwin T. Kruszewski 
ABSTRACT 

This dissertation deals with an analysis of impulsively loaded 
plastic plates that includes effects of large deformations and midplane 
forces . Specifically, it deals with a circular plate of uniform 
thickness simply supported at its edges. The impulsive loading is 
characterized by an initial velocity distribution. 

The analysis assumes that the plate is an isotropic, rigid, ideally 
plastic material. Shear deformations and longitudinal inertia are 
neglected. Both bending and midi Li ne forces are considered. The 
midplane forces can be either applied prior to the impulse or generated 
by deformations of the midplane. Deformations of the cross section are 
based on the assumption that a line initially normal to the midplane 
remains straight and normal after deformation. Strains are nonlinear 
with respect to transverse displacements but contain only linear radial 
displacement terms. Finally, the Tresca yield criteria are used. 

An interaction equation is derived between the plastic moment and 
normal force. Hie relationship indicates that the greater the midplane 
force the smaller the required bending moment for plastic flow. When the 



midplane force reaches a maximum value, the plate no longer carries a 
moment . 

Equilibrium equations are derived for the motion of the plate which 
are consistent with the von Rarman plate theory. Governing equations 
are then developed for the three possible phases of motion. The initial 
Phase 1 includes a bending hinge that travels from the support to the 
center of the plate. Phase 2, which is initiated when the hinge reaches 
the center, continues until either the motion ceases or a portion of the 
plate becomes a membrane. Phase 3 described the motion of the membrane 
hinge from its initial point of origin to its final stopping point. 

These governing equations are solved for two types of bending- 
moment— midplane-force interaction relationships. One relationship is 
based on displacement of the neutral surface from the midplane surface. 
The second is based on the magnitude of the midplane forces. Both types 
of plastic plate behavior are examined. 

Plots of final central deformation for various applied midplane 
forces and impulses are presented that clearly illustrate that even 
small amounts of midplane forces have a significant effect on the final 
deformation of plastic plates . Numerical results of the large def orma- 
tional analyses are also presented. Plots are given showing the 
influence of the magnitude of the Impulse on bending hinge velocity, 
initiation of membrane behavior, location of the membrane hinge, and 
the final deformation of the plate. 

A comparison is made between results of the large deformational 
analysi s and experimental data . It shows excellent agreement . For the 



small and intermediate range of impulses the agreement is within 
experimental scatter. For the large impulses the calculated deformations 
are slightly conservative. Reasons for this deviation are discussed. 
Finally a critical examination is made of the various other suggested 
approaches that have possible application to the behavior of plastic 
plates . 
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V. INTRODUCTION 


The problem of the impulsively loaded plastic structure is of 
interest to designers in several diversified areas. As early as 1951 > 
(see ref. 1) the U.S. Navy undertook a study of this problem in order 
to provide information to its ship designers about proper torpedo 
protection. The purpose of this study vas to understand the response 
and to predict the damage done to steel plates of ships and submarines 
by underwater explosions. The development of high energy explosive 
devices that can be detonated above ground has caused the designers of 
buildings to become interested in the problem. Here the impulsive loads 
are generated by shock waves traveling through the atmc 'phere rather 
than water. Rather recently, solutions to this problem have been sought 
by the designers of spacecraft. 

The outer shell of a spacecraft must be designed to withstand the 
rather harsh environment of space. Meteoroids are one of the more severe 
environmental hazards. Meteoroids are extra terrestrial particles, in 
orbit about the sun, '.hi'-h travel at velocities as high as 210,000 feet 
per second. Because of their extremely high velocity, even the smaller 
meteoroids can do a considerable amount of damage. For example, the 
large hole in the 3/l6-inch aluminum plate, shown in figure 1, was made 
by the small l/l6-inch-diameter projectile (also shown in fig. l) at an 
impact velocity of only 30,000 feet per second. At the highest meteoroid 
velocity this small projectile could penetrate a plate over an inch and 


a quarter thick. 



p 



Figure 1.- Damage to a 3/l6-inch aluminum target by 
l/l6-inch-diameter projectile. 
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In order to alleviate thi6 damage, spacecraft designers use a so • 
called "bumper” concept. This concept is based on the fact that vher 
a high velocity particle impacts a thin sheet of material both the 
projectile and the sheet material are shattered and vaporized into a 
cloud of debris. The use of this exterior wall or bumper thus changes 
the design consideration of the main wall from a penetration problem 
to one dealing with the momentum exchange between this cloud of debris 
and the main wall. 

Until recently, two basic approaches were used for the analyses of 
impulsive loaded plates: a pure membrane approach (refs. 1 and 2) and 

a pure bending approach (refs. 5 and ^) . Neither gives an entirely 
satisfactory answer. The former results in calculated deformations 
that are less than those experienced by the actual plate, hence uncon- 
servative. The latter gives deformations that are much greater than 
actual and hence too conservative. 

Recently, N. Jones (ref. 5) linked the two types of behavior. 

He permitted the plate to carry both maximum plastic moment and maximum 
plastic normal force during the first stage of deformation and full 
plastic normal force in the last stage . This, therefore, results in 
final deformations that are less than those of the membrane approaches 
and hence should be in poorer agreement with experiment. 

The purpose of this investigation is to develop a method of analysis 
that includes both bending and midplane forces and permits both flexural 
and membrane behavior of the plate. The analysis will deal specifically 
with a simply supported circular plate with an initial uniform velocity. 


1 + 


The equilibrium equations will be those appropriate to large deflections 
of plates as discussed in reference 6. The yield condition used in the 
analysis has been used by many investigators for moment axial -force 
interaction (see, for example, refs. 7 and 8 ). This yield condition 
defines the moment midplane -force interaction by examining the stress 
distribution in the cross section. It permits both to occur simultan- 
eously; however, both forces can never be maximum at the same time. 
Furthermore, one cannot exist if the other is at its maximum. 


VI. REVIEW OF PERTINENT LITERATURE 


General Considerations and Background 
There exists in literature only a limited number of theoretical 
analyses dealing with impulsively loaded plastic plates. The majority 
of the work that does exist is based on the assumption that the material 
is an isotropic, rigid plastic material whose behavior is consistent 
with the Tresca yield condition. The impulsive loading is characterized 
in these studies by the assumption that the structure is subjected only 
to an initial velocity distribution. 

The first analytical study of a plastic plate with such an initial 
condition was made by G. E. Hudson (see ref. l). In this study, the 
author treated a simply supported plate as a membrane, neglecting all 
bending loads and bending displacements. Initially, the entire plate 
was assumed to be plastic and traveling at a uniform velocity of V Q . 

A hinge developed at the boundary and traveled normal to it. 

The assumed deformation of the membrane is shown in figure 2. The 
hinge, at a distance p from the center, separates the circular 
membrane into two regions : a flat circular inner region that moves 

with a uniform transverse velocity of V Q ; and an outer region that is 
assumed to be rigid and conical in shape. Thus, plastic flow takes 
place only in the central or inner region. The hinge, with its 
associated discontinuity of slope, then travels with a uniform velocity 
until the inner circle is reduced to a point. Once the hinge reaches 
the center, it is assumed to remain there until the plate ceases to move. 
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An evaluation of this pure membrane approach is found in reference 2. 
Using an analysis similar to Hudson, Frederick compared his calculated 
results with experimental data on steel plates. Some of these results 
are reproduced in figure 3. In this figure experimental and theoretical 
deflection profiles are compared. Note that the calculated deflections 
are somewhat lower in value than those found experimentally. This, of 
course, is due in part to the assumption of a rigid outer region and the 
emission of all bending deformations. 

Frederick also examined two other factors that could influence the 
theoretically calculated deformations of the impulsively loaded membrane: 
the inclusion in the analysis of work hardening effects and the con- 
sideration of the Von Mises yield condition in place of the Tresca 
condition. Neither factor, however, was found to be of great importance. 
As shown in figure 3> work hardening has little effect on the final 
deformation. In fact, it tends to reduce the calculated deflections 
and hence, for this analysis, causes a poorer comparison with experiment. 

In direct contrast to the approaches discussed in references 1 
and 2, which considered only the membrane response of the plate, 

A. J. Wang (ref. 3) proposed an analysis of the deflections of 
a plastic plate under impulse loading that considered only the small 
deflection bending response of the plate. His analysis, like those of 
references 1 and 2, was based on the Tresca yield condition. His 
resulting flow rules, however, were different. 

The plastic flow rules, used by Wang, were first discussed by 
Prqger in reference 9« Hopkins and Prager (ref. 10) applied 
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Ideal plastic: 



Figure 3*- Membrane deflection profiles 
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these rules to dynamic problems, and developed a criteria for the 
existence and subsequent behavior of discontinuities across a hinge. 

In this work the author pointed out that if both the plate 
displacement and velocity are continuous, the following conditions 
across the hinge must be maintained. For a stationary hinge circle, 
the slope across a hinge in the radial direction is discontinuous, 
as, for example, in the case of the statically loaded plate in 
reference 11. For a moving hinge circle, the slope in the radial 
direction must be continuous. The accelerations for this case, however, 
can be discontinuous. 

Consequently, the pure bending solution presented by Wang (ref. 3 ) 
differs from the pure membrane approach in references 1 and 2 in several 
respects. The slope in the radial direction is continuous across the 
hinge. Plastic flow occurs in both the inner and outer regions of the 
plate. Finally, the velocity of the hinge is not uniform. 

Through the use of an expensive experimental program, 

A. L. Florence (ref. 4) made a critical evaluation of the small 
deflection bending analysis. In this experimental program, 22 aluminum 
and 20 steel pistes were subjected to impulsive loadings. The final 
measured deformations of these plates were then compared to theoretical 
calculations based on an analysis of Wang. 

Some of the results of reference 4 have been reproduced in 
figures 4 and 5 . In figure 4 theoretical and experimental nondixnen- 
sionalized deflections are plotted as a function of a nondimens ional 
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radial position, £. These results show that, when only small deflection 
bending loads and curvatures are considered, the theoretical results 
considerably overestimate the experimentally obtained deflections. 

Results in figure 4 are for one value of impulsive loading. How 
these comparisons are affected by changes in the magnitude of the 
impulsive loading is shown in figure 5. Here both theoretical and 
experimental deflections at the center of the plate are plotted as a 
function of the magnitude of the impulse. These results show that the 
correlation rapidly deteriorates with increasing deflection. This is 
primarily due to the increased importance of membrane forces and 
midplane distortion. Florence concluded from these results that a 
full treatment of the problem requires inclusion of membrane forces. 

Overeat imat Ions of calculated deflections by a bending analysis 
were also found by Gerald May and Kurt Gerstle (ref. 12) for the 
case of statically loaded plates. Here, a finite element approach was 
used in which elements were considered to be either entirely elastic 
or plastic . 

Jones (ref. 5) realizing the importance of the interaction of 
bending and midplane forces, attempted to link the bending and membrane 
behavior of plates in an analysis of impulsively loaded, simply supported 
plates. Similarly to the other investigators, he divided the motion of 
the plate into two phases. 

The first phase consisted of the initial deformations in which the 
hinge circle traveled from the support to the center of the plate. The 

analysis of this phase assumed that the central portion could resist 
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both a plastic moment, Mp, and a plastic midplane force, Np. The 
magnitude cf each of the forces, however, was assumed to be large 
enough to cause plastic behavior by themselves. In the outer region of 
this Phase 1 motion, the midplane force was assumed to be the maximum 
plastic value while the plastic moment ranged ^rom full maximum value 
at the hinge to zero at the support. Jones also assumed that the motion 
of the hinge was the same as that obtained by Wang for the pure bending 
assumption . 

In Phase 2 of the motion Jones assumed that only membrane forces 
were acting, that is, Mp = 0. Consequently, in Jones’ analysis, the 
plate always carried the maximum plastic midplane force, Np. In 
addition for a large portion of its motion, it also carries a plastic 
moment , Mp . 

The question of plastic midplane force and bending -moment inter- 
action has been quite extensively examined for beams (see, for example, 
refs. 7, 13, and 14). (Reference 1*1- is an excellent survey of analyses 
for plastic behavior; it is restricted, however, to work done on beams.) 
This interaction has also been examined for statically loaded plates 
(see, for example, refs. 15, 16, and 17). The yield condition that 
has evolved from all these studies relates the plastic loads present 
at any point in a structure by the following relationship: 

m p h \ 2 

=- + I "f- ) =1 where Mq and N c are the maximum plastic moment and 

Mo \^o/ 

axial force. It is this yield condition that will be used herein. 
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Review of Literature In Associated Areas 

An excellent review of all the phases of plasticity essential to 
development of plasticity analyses is given by Olszak, Mroz, and Perzyna 
in reference 15. This summary deals not only with the fundamental 
basic approaches and constitutive equations but also stresses the 
application of the theory to various practical problem areas such as 
plates and shells, soil mechanics, three-dimensional problems, and 
axially symmetric problems . This reference also devoted a whole chapter 
to dynamical problems . Efech section is followed by a complete 
bibliography. 

An equally thorough survey, but limited only to analyses of beams 
under dynamic loading, is found in reference 1J. 

The preceding section dealt only with the state of the art 
pertaining to the specific problem being analyzed herein; that is, a 
simply supported circular plate constructed of a perfect rigid plastic 
material and subjected to an initial uniform velocity. There are, 
however, analyses dealing with other problems that need mentioning 
both for the sake of completeness and, in some cases, for their possible 
application to the large deformation aspect of the present analysis. 

In references l8 and 19, Thomson studied the plastic behavior of 
simply supported circular plates under nonuniform transverse impulses. 
Only the bending behavior of the plates was considered. The plate was 
assumed to be impacted by an axisymmetric impulse that induced 
an initial velocity with a general radial Gaussian distribution. 
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Perzyna (ref. 2o) investigated the effects of a time variation 
of a loading on plastic deformations of a simply supported circular 
plate. This aralysis, like Thomson's, assumed that the material was 
a rigid plastic behaving in a manner consistent with the Tresca yield 
criteria, and considered only the bending mode of deformations. The 
pressure loading was taken to be uniform in both the radial and 
tangential directions. The magnitude of the pressure, however, varied 
with time. Perzyna found that for a given total impulse the character 
of the time pressure function had little effect on the final shape of 
the plate. 

Plastic behavior of plates with other types of boundary conditions 
has also been investigated for the bending case. Wang and Hopkins 
(ref. 21) gave results for circular plates vith built-in edges. Shapiro 
(ref. 22) investigated the plastic response of an annular plate. This 
rigid plastic plate was assumed fixed at inner radius and the outer rim 
was subjected to a constant velocity for a finite length of time. 

Finally Mroz (ref. 22) investigated the case of simply supported 
annular plates for two different conditions : for a uniform lateral 

loading with simply supported inner ring and for a uniform shear force 
applied to the inner periphery whose edges are considered free. 


VII. SIMBOIS 
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F(n) 

I 
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^0 
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Qr»% 

R 

Vo 
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nondimensional central deformation of the plate 
value of A when plate motion stops 
constants of integration (see eqs . (103), (H5), and 
(129)) 

function of n defined by equation (131) 

nondimensional impulse defined in equations ( 16 ) and (17) 

Bessel function of the first kind of order 0 

curvature of plate in r and 8 directions, respectively 

bending-moment resultant 

value of M for plastic behavior 

value of M at the hinge circle 

2 

mnYlimim bending moment plate can sustain, o 0 h 
value of M in the r and 0 directions 

twisting moment 
midplane force 

value of N for plastic behavior 

value of N at the hinge circle 

maximum midplane force plate can sustain, 2a 0 h 

value of N in the r and 0 direction 

shear force in r and 0 direction 

radius of plate 

initial velocity of plate 

nondimensional displacement of plate, ^ 
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17 


function defined by equation (46) 

nondimensional distance between midplane and neutral 

_ z 0 
surface. — 
h 

half depth of plate 

natural logarithm 

nondimensional membrane radius 

initial value of n 

final value of n 

loading per unit area of plate 

radial coordinate 

distance along midplane in radial direction 
dimensional time 

radial displacement of middle surface 
transverse deflection of middle surface of plate 
coordinate along axis of symmetry 

distance between middle surface and neutral surface of 
the plate 

angle between the neutral surface and r coordinate 
strains at any point of cross section in r and 9 
directions, respectively 

strains in midplane surface in r and 0 directions, 
respectively 
dummy variable 
dummy variable 


tangential coordinate 
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X 

M 

5 

p 

P 1 

a r> a Qt a z 

a o 

a I* a II* a III 

T rd’ T z9 ,T rz 

T 

T 1 

t 2 

Subscripts 

r 


constent 

mass of plate per unit area 

nondimens lonal radial coordinate., ~ 

R 

plastic bending hinge circle radius 
value of p at the initiation of the membrane hinge 
normal stress in r, 9, and z direction, respectively 
yield stress 

principal stresses (Jj > Cjj > a-j-jj 
shear stresses 

V 0 t 

nondimensional time, 

value of t when plastic bending hinge leaches center 
of plate 

value of t vhen plastic membrane hinge is initiated 
radial direction 


0 


circumferential direction 


t differentiation with respect to t 

t differentiation with respect to t 

A line over the symbol M denotes a dimensional moment such that 
^ 

M = — . A line over the symbol N denotes a dimensional midplane such 

— n 

that N = — . A dot over the symbol denotes differentiation with 
«o 

respect to time , t . 


VIII. DERIVATION OF BASIC EQUATIONS 


Assumptions 

The purpose of the present work Is to develop an analytical method 
for the prediction of the plastic behavior of impulsively loaded plates 
that will combine both membrane forces and bending. To do this the 
analysis must contain large deflection effects and permit midplane 
distortions along with bending distortions. 

The plate that is to be analyzed and the coordinate system that is 
to be used is shown in figure 6. The plate is circular, has a radius 
R, a thickness 2h, and is simply supported at its edges. The imposed 
initial conditions are that the displacement, w(r,0), is zero and the 
initial velocity, w-t(r,o), is a constant, V Q . 

The analytical method will be based on the following assumptions. 

The plate is assumed to be an isotropic, rigid, ideally plastic material. 
Consequently, both elastic and work hardening effects will not be 
included in the analysis. Although midplane forces and distortions 
are permitted, inertia effects in the radial direction are omitted. 

Shear deformations and rotary inertia effects are also neglected. The 
equilibrium and resulting governing equations used in this analysis will 
be consistent with those of the von Karman plate theory. This theory, 
as discussed in reference 24, neglects third and higher order terms of 
the dependent variable so that 


sin 0 *» 3 
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and 

cos 3 *= 1 - § 3 2 

The following sections will discuss some of the basic governing equations 
for such behavior. 


Equilibrium Considerations 

The forces acting on the plate and on a typical element are shown 
in figure 7. These forces consist of bending moments Mj. and Mq, 
midplane forces N r and N 0 , a shear force Qj. , and a uniformly 
distributed loading p. In the case of the loading symbols, all symbols 
with bars above will denote dimensional forces, while the force symbols 
without the bars will represent nondimens ional forces. All forces are 
shown in figure 7 in their positive direction and all are specified per 
unit length of line element in the middle surface of the plate. Due to 
the cylindrical symmetry of loading, the shear stresses T r0 and t z0 
vanish . Consequently the shear forces Qg and the torque Mg r also 
vanish. With the assumption that the plate thickness (2h) is much 
smaller than the radius R, the influence of stresses a z and T rz 
is small and hence negligible. 

From figure 7(a) the summation of forces in the midplane direction 
results in 




= 0 


( 1 ) 
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From figure 7(b) the summation of forces in the transverse direction 
gives an expression for the transverse shear, Qr, 

Or = -N r ^ j (P - UW tt )C dC (2) 

where the term uw tt represents the inertial loading due to transverse 
motion, and u is the mass of the plate per unit area. The problem 
discussed herein deals only with the impulsive loading that results in 
an initial velocity of the plastic plate. The loading term, p, 
therefore, is zero and will not be inclu^jd in the remainder of the 
analyses . 

The summation of moments in figure 7(a) results in the equation 

dM r _ _ 

Mj. + r ~ - Mq = rQj. (3) 


Eliminating the shear force Qp from equation ( 3 ) (see eq. ( 2 )) 
results in the equilibrium equation 


|;(rM r ) - Mq = 




n^ttC 


(*0 


Equations (1) and (k), therefore, represent the equilibrium 
equations for the plastic plate with midplane forces. These equations, 
of course, are identical to the vo.i Karman plate equations obtained by 
considering large deflections of circular plates (see, for example, 
ref. 6) . 


Nondlmenslonalization 


In order to simplify the derivations and to obtain the natural 
parameters of the mechanism involved, the remaining analyses will be 
written in nondimens ional terms. The following nondimens ional parameters 
are therefore defined. 

The moment is normalized with respect to Mq where 

M 0 = h 2 a 0 (5) 

and <j 0 is the critical yield stress. More will be said about this 
stress in the section on Plastic Flow Considerations. Mq represents 
the maximum bending moment the plate can possibly carry. Consequently, 
the absolute value of the nondimens ional i zed bending moment 


M 



(6) 


can never be greater than unity. Similarly the nondimensionalized 
normal force is 


N 



(7) 


where 
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N 0 is the maximum normal force possible and, similarly to M, | W | 
can never exceed unity. The nondiaensional form of the shear force, 

Qr, is 


Or 



(9) 


The nondimensionalized radial coordinate, f, displacement, W, 
and time, t, are defined as follows 



( 10 ) 


w 



( 11 ) 


and 



( 12 ) 


Equilibrium equations (1), (2), and (4) can now be written in 
nondimens ional form as 


^(|N r ) - N 0 = 0 


(15) 


Or = 


-N„ 


+ 51 



(14) 
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J~( tMj.) - Me - -2|B r * 61 J* W TT * d» 


where I is a nondimens ionai ized momentum defined as 


uV 0 2 R 2 

12h 2 N„ 


2 2 


Deformation Considerations 

Consistent with all of the previously mentioned analyses on the 
plastic behavior of plates, the deformations of the cross section will 
be based on the assumption that a line initially normal to the midplane 
remains straight and normal after deformation* Consequently the strain 
at any point in the cross section of the plate can be written as : 

For the radial direction 


e r = e r + zK r 


For the tangential or 9 direction 


€ e = € e + zK e 
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where u is the displacement of the midplane in the r direction 
(u is considered positive when the displacements are in che positive 
r direction). If only small displacements are to be considered, the 
second term in equation (22) can be neglected. Note that the strains 
are written from a Lagrangian viewpoint and no nonlinearities due to 
radial displacements are considered. 
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Plastic Flov Considerations 

Plastic behavior of the plate will be governed by the Tresca yield 
criteria. This criteria states that plastic yielding will begin when 
the maximum shear stress reaches a certain critical value, that is 


= q I ~ q III = fo 
T max 2 2 


(24) 


where <jj and Ojjj are, respectively, the maximum and minimum 
principal stresses and a Q the critical yield stress of the material. 
With the assumption that the transverse normal stress is negligible 
(c z = 0) , the Tresca condition can be graphically represented as shown 
in figure 8 in terms of a T and cr 0 . 

Any combination of dp and cr 0 that falls on this boundary will 
initiate plastic flow. Any subsequent plastic flow will occur tangent 
to the surface of the yield diagram. For example, if cr r and Og 
fall on the line AB (see fig. 8), plastic flow will occur perpendicular 
to AB in a manner that the radial strain rate, l r , is zero. For 
stress combinations along AF, 1 Q will be zero. At the vertices of 
the diagrams the specific direction of plastic flow is not defined; 
both stresses, however, are defined. Thus, for all points on this 
diagram two conditions are specified: either both stresses or a 

direction of strain rate and a relationship between the stresses. 

For pure bending of a plate the moment along a plastic flow line 

p 

is unique and equal to h 0 O . Similarly when only midplane forces are 
considered the plastic membrane force is unique and equal to 2h a Q . 
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However, when both midplane forces and bending are permitted, then 
there are no unique values of M or N at the initiation of plastic 
flow. 

A typical cross section for this case is shown in figure 9* 
Notice that the development of a midplane force causes the plastic 
stress distribution to create a neutral surface at some distance z Q 
from the midplane. The plastic moment and membrane force are then 
given by 


M p = (h 2 - z 0 2 )a 0 

and 


(25) 


w p = - 2z o a c 


(26) 


Notice that the plastic moment and membrane forces are coupled through 
z o • 

This coupling indicates that, while there is no unique value of the 
plastic midplane force, N p or moment M p , there is an interaction of 
the two. Elimination of z Q from equations (25, and (26) gives the 
form of this interaction 


Mo 



= 1 


( 27 ) 


In nondimens ional form equations (25), (26), and (27) can be 


written as 
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Figure 9.- Distribution of stresses in cross section of plate. 
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1 - Z 


N p = -Z 


M p + V = 1 


z= f° 

Z h 


The interaction equation (eq. (27) or eq . (30)) gives a relationship 
between the plastic moment Mp and the plastic normal force Np. As 
long as Np is less than N 0 (Np < 1) a plastic moment exists. However, 
as soon as the plastic normal force becomes N Q (Np = 1) then Mp 
becomes zero and the section no longer carries moment. These interaction 
relationships are especially useful if the value of the lateral load is 
prescribed, as in an analyses involving the bending behavior of a 
plastic plate with applied lateral loads . 

When the deformations of the plate become large, lateral loads are 
generated by the stretching, u(r,t), of the midplane. 

With the condition that at the support u(R,t) is zero, an 
integral expression for u can be written 


(ds - dr) 
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where s is the distance measured along the midplane, 
sketch it can be seen that 


From the above 


ds 2 = dw 2 + dr 2 


(33) 


Substitution of equation (33) in equation (32) results in 


u = 



which, for the case where 


dw 

dr 


< 1 , 


can be approximated by 


( 3*0 



(35) 


The separation of the neutral and midplane surfaces creates a 
coupling between slope and midplane distortion, u. This coupling is 
achieved through the previously mentioned assumption of straight normals . 
That is, a line perpendicular to the neutral surface before deformation 
will remain straight and normal after deformation. Consequently, in the 
radial direction the midplane distortion, u, is related to the radial 
slope by the equation 


With equations (35) and (3 6 ) z 0 can be expressed in terms of the 


slope as 



dr 


or in nondimens ional terms 

5 £ (I h m 

z = r 

dw 

*1 

The magnitude of Z depends on the shape of the displacement . 

From equation (29) and the knowledge that j Np J < 1, |z| must also 
be less than 1. Consequently, those portions of the plate where the 
absolute value of Z is greater than unity must be acting as a membrane 
(that is, they carry no moment) and there is no valid unique relationship 
between slope and displacement. 


IX. DERIVATION OF GOVERNING EQUATIONS 


Discussion of Physical Behavior 

The following section Is a description of the physical behavior of 
a plate as it reacts to the initial Impulse. It should not in any way 
be construed as a statement of assumptions. On the contrary, the 
mechanisms and behavior described all result from the basic assumptions 
and considerations discussed in the previous chapter. These will become 
evident as the analysis in the following sections proceeds. 

The purpose of this approach is twofold: First, it will give the 

reader a clearer understanding of the motivation and significance of 
the analytical work that follows. Second, it will permit the author to 
define and classify certain regions and behaviors before proving their 
existence. 

Just after the initial moment of impact, the plate is separated 

into two plastic regions by a hinge circle, as shown in figure 10 . 

The hinge, which originates at the support, travels at a velocity 

P T until it reaches the center. The central portion of this plate 

(0 < r < pR) travels at a uniform transverse velocity equal to the 

initial velocity, V Q . Such motion has been observed experimentally 

(see ref. 26). The stress conditions are such that a = <j„ = a . At 

r 0 o 

any time, t, the displacement w in the central region is V Q t . 

In nondimens ional terms the central displacement is 

W(0,t ) = A(t) = 2t 0 < 5 < p (39) 
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Thus the nondimensional velocity in the central region is 

W t (0,t) = A^t) =2 0 < { < p (4o) 

In the outer region (p < | < 1) only the tangential stress, Oq, 
is plastic. The velocities are linear w.'th 

W T (p,r) = 2 (U) 

and 

W t (1,t) = 0 (k2) 

This phase of the motion will be referred to as Phase 1 . Note that 
the displacement of any point on the plate during this phase depends 
only on the function p(r) . Thus the parameter p could be thought 
of as a time parameter. 

Once the hinge circle reaches the center of the plate. Phase 2 of 
the motion begins (see fig. 10). This phase of motion is similar to 
that occurring in the outer region for Phase 1. The Gtress, Oq, is 
plastic and the velocity distribution linear. 

If during either Phase 1 or Phase 2 the motion and loading of the 
plate is such that midplane force N becomes 1, the behavior of the 
plate is changed. 

A hinge circle, which will be referred to as a merbrane hinge, 
forms at thi 6 point (5 = n) (see fig. 10). This hinge divides the 
plate into two regions. The outer region, n < | <1, carries both 
midplane forces and moments and behaves as before. The inner portion 



10.- Three possible phases of motion of plastic plates. 
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0 < | < n behaves like a membrane. It carries no bending and transmits 
only a constant tensile load. This membrane hinge then travels at a 
velocity, n y , toward the support. This phase of the motion will be 
referred to as Phase 3* 

Analysis of Phase 1 Behavior 

This phase of the plate behavior is typical of the initial phase 
of motion used by all of the previously mentioned investigators of 
ideally plastic, circular plate behavior (see, for example, refs. 10 
and 11). Because of the rotational symmetry of the loading the maximum 
stress in a circular plate must occur at the center. Thus the plate 
cannot reach the flow limit without becoming plastic at the center. 

Due to rotational symmetry the radial and tangential stresses must be 
equal, thus a r = Og = cj 0 . Therefore, the material at this point is in 
regime A of the T* 'sea yield hexagon (see fig. 8). At the simply 
supported boundary, the radial moment must vanish so that a r = 0. 

This portion of the plate must be in regime AB of the Tresca yield 
hexagon . 

Consequently, the plate is assumed divided into two regions : & 

central portion in which the material is in regime A, and an outer 
region where the material is in regime AB. Biese two region? are 
separated at | = p by a plastic hinge. 

The inner portion of the plate < p) forms a circular plate with 
midplane forces at its outer edge (£ = p). The magnitude of these radial 
midplane forces, B r , is .aken to be N p . Because of the symmetrical 
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loading condition imposed on the plate, the stresses at the center must 
be such that N r = N • A solution that satisfies these boundary condi- 
tions and the equilibrium equation (eq. (13)) i6 that the stresses are 
such that N r = Nq for all interior points. Therefore, N is not a 
function of position and 

Np = N e = N p 0 < | < p (43) 

for all points in the inner region. Since the midplane force, N p , 
can be a function only of time, the plastic moments must also be a 
function only of time (see eq. (30) so that 

M r =M 0 =M p 0 < | < p (44) 

With equations (43) and (44) the equilibrium equation (eq. (4)) reduces 
to 


^ r I 

6l / w t dC = 2|N — o < I < P '45) 

J o p d 5 

In order to solve this governing equation a separation of variables 
technique will be used in which 

W(|>t) = A(t)w(| ) (46) 

Substitution of equation (46) into the integral differential equation 
(eq. (45)) results, after some simplification, in 


4o 


du 

H 


J i d£ 


31^1 

dr 2 


N p A 


(47) 


Since the left side of equation (47) is a function of | and the right 
side a function of t, both must be equal to some constant. A, 


dw 

1 a* 
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a 2 A 


f 5 - 


V 

\pT 


= -A c 


W5 


(48) 


Hence the governing equation for W is 


| ^ + A 2 f 1 W£ d£ = 0 
°l J o 


(49) 


Differentiating equation (49) gives 


, d c w . dw . ,2rr t 

*7? * 1 


= 0 


(50) 


Solution of equation (50 ) results in 


W = J 0 (A|) 


(5D 


The deflection in this inner region can now be written as 


W(|t) = AJ C (A|) 


0 < P < | 


(52) 
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The initial condition of the problem states that the velocity 
everywhere is equal to 2 so that 

w T (s,o) = A t j c (M) = 2 (53) 

Thus A must be zero and the velocity in the inner region is 

W T = 2 0 < | < p (54) 

so that 

W TT = 0 o < I < P (55) 

W = 2t 0 < | < p (56) 

In the outer region (p < 5 < l) the following conditions apply. 

At the hinge, | = p, the stresses are continuous, so that (%. = Oq = 0 o . 
At the support, g = 1, 0j> = 0. This portion of the plate, therefore, 
is in regime AB of the Tresca yield criteria. In this regime the 
plastic flow is such that the strain rate in the radial direction T r , 
is zero for all points in the cross section (see eq. (l8)) so that 

? r = e r + zK r =0 (57) 

Since this equation mu6t hold for all values of z 


€ r = K r = 0 


(58) 



This condition coupled with equation (20) gives 


d^W 


= 0 


( 59 ) 


As the velocity is continuous across the hinge 

W T (p,r) = 2 (6o) 

While at the support | = 1 

W T (l,x) = 0 (61) 

The solution of equation (59) consistent with boundary conditions, 
equations (6o) and (6l), gives 


"t - 2 (ri^) p < t < 1 < 62 > 

Equation (62) represents the velocity distribution of the plate in 
the outer region. Since p is a function of t the acceleration of 
the plate in this region is 



= 2 


(1 - I) 

(1 - P) 2 


P < 5 < 1 


( 65 ) 


where p^ is the velocity of the hinge circle. 

A relationship for p T can be obtained through the use of the 


equilibrium expressions. 
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Integration of the midplane-force equilibrium equation (eq. (13)) 
with the boundary condition that N r = Nq = N p at the hinge (see 
eq. (43)) gives 


|N r = 



(64) 


Integration of the moment equation (eq. (15)) with the boundary 
conditions for moment at the hinge, Mp = Mq = M p , (see eq. (44)) 
gives 


|Mr = pM p + 



^ d£ - 2 







d \| f 


(65) 


Ifce third term on the right side of equation (65) can be rewritten 
through the use of integration by parts as follows 



d| = |WNj. 


r 


w 


6(SN r ) 




(66) 


This integral can now be written in terms of N p and Nq through the 
use of equations ( 13 ) and (64) 



pWN c 


V /" 5 N 0 dt, - f WH 0 d| 


(67) 


Substitution of equation ( 67 ) in equation (65) gives 


|Mj. = pM p + f 1 Mq d£ - 2pWN p + 2pW(p)N p - 2W [ 5 N 0 d£ 


+ 2 r 6 wN e d? + 61 [* u rr t 

v * J p o 






(68) 


Equation (68) gives the value of the radial moment at any point 
The governing equation for p T is obtained by evaluation of equation (68) 
at the support | = 1. 

At the support 


W(l,r) = 0 (69) 

and 

Mj.(1,t) = 0 (70) 


Thus 


P M n + 2pW(p)K- + 






d\jr 


(71) 


Through the use of equations ( 55 ) and (65) the inertia term of 


equation ( 71 ) can be evaluated 


*+5 


Substitution of equation (72) into equation (71) gives 


1(1 - p)(l + 3 p)p t = -pM 


2pW(p)N p - 


Me ae 


- 2 


WM 0 d 5 


(73) 


From the interaction equation between plastic normal forces and 
plastic moments (eq. (30)) and the expression for W(p) from 
equation (56) (W(p) = 2t), equation (73) becomes 

1(1 - p)(l + 3 p)p t = -1 - 4pN p r + pN 2 + I N 0 2 d| - 2 [ VN 0 d| 

J P J P 


(7*0 

Equation (7*0 can also be written in terms of the displacement of 
the neutral surface with the use of equation (29) 


1(1 - p)(l + 3 p)p t 


•1 + ^pZt + pZ p 2 + 




WZ d| 


(75) 

Equations (7*0 and (75) are the governing equations for the hinge 
circle in the Riase 1 motion. The equations are dependent on the shape 
of the displacements and the values of the midplane forces. If only 
bending is considered, that is Nq * 0 or Z * 0, the equation reduces 
to that obtained by Wang in his bending solution for impulsively loaded 
plates (see ref. 3). 
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Analysis of Phase 2 Behavior 

Phase 1 motion will continue until one of two conditions occurs . 
Either the hinge circle reaches the center of the plate or a portion of 
the plate becomes a membrane (that is, the plastic moment vanishes). 

The former case will be considered in this section. 

Let be the time it takes the hinge circle in Phase 1 to reach 

the center of the plate. Then the displacement and velocity at that 

time will be (see eqs. (5*0 and (56)) 

W(0, T;l ) = 2 t x (76) 

and 

V T (0, Tl ) = 2 (77) 

Hence these conditions form the initial conditions fpr-the Phase 2 
motion. 

The radial stresses in this phase range from at the center 

to zero at the support. Hence the material behaves in a manner similar 

to that in the outer region of the Phase 1 motion. The radial strain 

rates are zero and the velocity varies linearly from A T at the center 
to zero at the support. Thus 

W T * A T (l - |) (78) 

and 

W tt =A tt (1-!) (79) 

where A represents the deflection at the center of the plate. 
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The governing equations for this phase are identical to 
in Phase 1 with p set equal to zero. 


6l 


[ 1 I” [ * W TT g dg d* = - WN e c 

° o L J o -J ^ o o 

The inertia term for this phase becomes (see eq. (79) 

61 [ l / w g dg d* = I ia tt 
^o L' J o J 

Substituting equation (8l) into equation (8o) gives 




r 1 

Mq d| - 2 / WNq d| 


which in terms of Ng is 


| IArr - -1 + f 1 » e 2 d| - 2 f 1 V K e d| 
* J o d o 


and in terms of Z is 


h IA tt - -1 + f Z 2 d| + 2 f 1 WZ dg 
J o J o 


equation (71 ) 


(8o) 


(81) 


(82) 


(83) 


(84) 


Equations (83) and (84), as in Phase 1, reduce to the pure bending 
case for N Q = Z = 0 . These governing equations apply until either the 
motion of the plate stops or a portion of the plate becomes a membrane. 



Equations ( 7*0 and (or eqs. ( 75 ) and ( 84 )) are sufficient to 

define the motion of plates where the applied impulse, I, is small. 

When the impulse becomes large enough to cause a pc --t ion of the plate 
to be a membrane, an additional phase of motion must be considered. 

Analysis of Phase 3 Behavior 

As the deformations Increase the midplane forces also increase. 

If the initial impulse, I, is large enough, the midplane forces In the 
central portion of the plate will become equal to one. Thus thf. bending 
moment vanishes and this portion of the plate behaves as a membrane. 

Let the time at which a portion of the plate becomes a membrane be 
t 0 . As time increases, the size of the membrane region also increases. 
The plate is thus separated into two regions : one in which there is no 

moment and a second in which there i6 an interaction of moment and 
midplane force. 

In the inner region, 0 < | < r., Mq = 0 and N 0 = 1. In addition, 
N r = N 0 at the center. These conditions along with equilibrium equation 
(eq. (15)) result in 

N r = N Q = 1 0 < | < n (85) 

for the entire inner region. Consequently, this region is in regime A 
of the Tresca yield diagram. The point of separation of the membrane 
portion of the plate, 5 = n, can then be thought of as a membrane hinge. 
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The moment equilibrium equation (eq. (15)) reduces at | = n to 


n 


8w 

* 


r n 

= 31 / v TT | d* 

J o 


( 86 ) 


is evaluated at | = n. 


n 


dw 

where the slope r— 

The governing equations for the motion of this hinge can be obtained 
from the governing equation for the outer portion of the plate 
(n < | < 1). Hiis governing equation is identical to that for outer 
portion of Phase 1 motion with n substituted for p (see eq. (71)). 


2nW(n) + 




d ♦ 


( 87 ) 

(Note that the conditions N p = N n = 1 and M p = Mjj = 0 were vied.) 
With the use of equation (30), equation (87) becomes 




d* = -(1 - n) 


- 2nW(n) + 


r 

d n 


W0 2 d| - 2 



(88) 


In terms of Z (see eq. (29)) 




W 


TT 



di|r = -(1 - n) 




n < C < 1 


( 89 ) 
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As in Phase 1 motion, W|| T = 0, so that 


■ <H T> n (rri) 


n < 5 


where (W ) is the velocity at the membrane 
T n 


the velocity is linear 
< 1 (90) 

hinge. 


X. SOLUTIONS OF GOVERNING EQUATIONS 


In Chapter IX the governing equations for the bending hinge, p, 
the membrane hinge, n, and the center deflection. A, were derived. 

No assumptions were made in this section so that these equations are 
exact within the framework of the basic assumptions, equilibrium 
equations, and flow conditions discussed in Chapter VIII. 

Note that these governing equations contain either a Ng or Z 
and a deflection W. Consequently, a solution can be obtained if seme 
assumptions concerning these terms are made. 

The deflection term appears only in a definite integral. Further- 
more, the value of the deflection is always known at both endpoints of 
the definite Integral, that is, W(p) * 2 , W(o) = A, and w(l) = 0. 
Consequently, only the assumption of shape of the deflection profile 
is necessary. 

An assumption of conical surfaces will be made in the sections to 
follow. Such deformat ional surfaces have been observed experimentally 
for dynamic loading (see ref. 26) and have been used by others in 
solving membrane type problems . For example, both Hudson in reference 1 
and Frederick in reference 2 made this assumption in their treatment of 
impulsively loaded membranes. Onat and Haythornthvaite (ref. 25) also 
made the same assumption in a study dealing with the influence of midplane 
forces on plastic plates under static loading. 

In order to solve the governing equations in the form of 
equations (7*0 and (83), some knowledge of the value of Ng is also 
necessary. If one recalls that in the elastic solution a circular plate 
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uniformly loaded at its supports has constant midplane stresses throughout 
the plate, it is tempting to assume that the plastic plate has a similar 
distribution. Applying this reasoning to the plastic plates with large 
impulses, the entire plate could be assumed a membrane and N 0 set 
equal to one. 

This approach, hovever, does not take into account midplane forces 
that are developed by the distortion of the midplane of the plate. The 
failure to recognize this fact reflects itself in two ways: First the 

midplane forces are a constant. Second, and more Important, no portion 
of the plate is permitted to become a membrane. Thus an accurate 
solution to impulsively loaded plates must be based on the governing 
equations that contain the neutral surface displacement Z (eqs. (75) 
and (84)). it also must, when appropriate, utilize the Phase 3 governing 
equations (eqs. ( 89 ) and (90)). Solutions using the assumption of 
constant midplane forces are of interest. They can serve as an 
evaluation of the influence of midplane forces on the deflections of 
plastic plates. In addition, the solution for N Q = 1 gives a limiting 
case of maximum midplane force. Thus both solutions will be developed 
in the following sections. 

Bending of Plates Under Lateral Loads 

In the analysis to follow the midplane forces will be assumed 
constant throughout the plate. With the exception of the limiting case 
of N q = 1, the plate will always carry a plastic moment and at no time 
will the plate become a membrane. Consequently, this case can be thought 
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of as an analysis of the bending behavior of a plastic plate that has 
been loaded by lateral loads, N Q , prior to the application of the 
impulse. This analysis will not recognize midplane deformations due 
to large deflections . 

The governing equation for Phase 1 (eq. (7*0) reduces for this 
case to 


1(1 - p)(l + 5 p)p t = N 0 2 - 1 - *+pN 0 r 



(9D 


If the deflected surface is assumed to be a conical surface 


W = 2t 




(92 ) 


so that w(p) = 2t and W(l) = 0. Substitution of equation (92) into 
equation (91) results in a nonlinear differential equation 

1(1 - p)(l + 3 p)p t = -(1 - N 0 2 ) - 2N q t( 1 + p) (93) 

This equation is solved in appendix A through the use of a finite 
difference approach. Numerical results are given for values of N 0 
ranging from 0 to 1.0 and for a range of Impulses, I, from 0.*+ to 11.6. 
Exact solutions are obtained for the limiting cases of Nq equal to 
zero and 1. Time histories of the hinge circle, p, are tabulated in 
appendix A for a widq range of I and N 0 . 


Phase 1 motion will continue until the hinge circle reaches the 
center. Let be the time at which this occurs. At t = the 

deformation at the center will be 

wCo,^) = ACtj,) = 2t 1 (9*0 

The velocity at the center will be 

W t ( 0 ,t 1 ) = A T ( Tl ) = 2 (95) 

The value of is tabulated in appendix A as a function of impulse 

I and midplane force N 0 . The values of t 1 with the corresponding 
values of displacement A(t^) and velocity 2 serve as the initial 
values of the Phase 2 motion. 

The governing equation for Phase 2 (eq. (85)) reduces to 

| IAtt - "(I " He 2 ) - 2% f 1 W d{ (96) 

J o 

Consistent with the assumption of conical surfaces 

W = A(1 - |) (97) 

Thus the governing equation becomes 

IAtt « - 2 U - N 0 2 ) - 2N 0 A (98) 

The first integration of equation (98) can be achieved by 
multiplying the equation by A T . 
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“At - -(?U - »e 2 ) * 2 »e*k 

(99) 

The integral of equation (99) can now be written as 


| IAt 2 - -2(1 - Nq 2 )A - HqA 2 + Ci 

( 100 ) 

where Ci is the constant of integration. 
The initial conditions for Phase 2 are 


A(t x ) = 2 Tl 

( 101 ) 

A t (t 1 ) = 2 

( 102 ) 

Therefore 


C-l - 21 + 4 Tl (l - N 0 2 ) + 4 Tl 2 N 0 

(105) 


Substitution of equation (105) into equation (100) gives the 
expression for the velocity of the center of the plate a 6 

HU 2 = -Ml - Nq 2 )A - 2N 0 A 2 + 4l + 8^(1 - N 0 2 ) + Bt^Nq (loM 

The final deformation at the center of the plate. Ay, can now be 
obtained from equation (lo4) by setting A^. = 0. 

N Q Ap 2 + 2(1 - N 0 2 )A f = 21 + 4 Tl (l - N q 2 ) + 


(105) 
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Large Deflection of Plates 

The analyses in this section will consider the midplane forces 
that are developed by large deformations of the plate. Consequently, 
governing equations (75 ), (84), and (89) will be used as these equations 
contain the parameter Z. As previously discussed, this parameter is a 
measure of the displacement of the neutral surface from the midplane. 

The definition of Z is given in equation (58) in terms of the 
midplane displacement W. For the assumption of conical surfaces 
equation (38) reduces to 


Z * 


W 

2 


(106) 


For Phase 1 


and 


Thus 


W = 2t 


0 < | < p 



P < t < 1 


Z = -T 


0 < I < p 




p < i < 1 


(107) 


(108) 
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The governing equation for Riase 1 can now be written as 

ia - p)(l + 3 p)p t - -1 - (2p + 1 )t 2 (109) 

A numerical solution of equation (109) using a finite difference 
approach is given in appendix B. Numerical values are calculated for a 
range of Impulses from I equal to 0.2 to 12. 

Note that at t = 0 equation (109) reduces! to the differential 
equation for pure bending (see ref. 3)* ®ius for small values of r 
the behavior of the hinge approximates that obtained for a plastic plate 
in bending. As time becomes large, however, the magnitude of the hinge 
velocity p T increases over that for the bending case. This increase 
in hinge velocity results in a decrease in final deformation. 

Phase 1 motion continues until the hinge circle reaches the center 
of the plate. At this time, t » t^, deformation of the center is 

V(0,t 1 ) » A(t 1 ) - 2 t ± (110) 

and its velocity is 

W T (0, Tl ) - A T ( Tl ) = 2 (111) 

Thus equations (110) and (ill) serve as initial condition.* for 
Phase 2. The governing equations for this phase are given by 
equation (84). Using equation ( 106 ), equation (84) becomes 

I I*TT * -1 * £ 


( 112 ) 
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which can be approximated with the use of the conical surface assumption 
by 

2IA tt * -4 - A 2 (115) 

By a procedure similar to that used for the solution of 
equation ( 98 ), equation ( 115 ) can be integrated to get 


IA ^ 2 = - |(12 + A 2 ) + C 2 (114) 

The constant of integration C 2 can be evaluated by use of the 
initial conditions (eqs. ( 110 ) and ( 111 )) 

C 2 * 41 + “T (l2 + Ut 1 2) (115) 

Substitution of equation (115) into equation (114) gives an 
expression for the velocity of the center of the plate during Phase 2 
motion 


IA t 2 = 4l + {l2 + (2r x ) 2 } - j(12 + A 2 ) (116) 

The final deformation of the center Ap can be evaluated from 
equation (116) by setting A T equal to zero. Thus 


A F (12 + Ap 2 ) - 121 ♦ (2t x ) (l2 «■ (2 Tl ) 2 } 


(117) 
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At first glance it would seem that equation (117) can provide a 
solution for the fined deformation of the center of all impulsively 
loaded plates. There are, however, restrictions on the use of this 
equation. 

Equation (117) is based on the assumption that only Phase 1 and 
Phase 2 behavior occurs during the deformation of the plate. Both of 
these phases assume that an interaction of bending and midplane forces 
occurs. It has been shown that as the plastic midplane forces, N, 
increased the value of the plastic moments, M, decreased until a 
value of N equal to one is reached. The moment at that time is zero 
and that portion of the plate became a membrane . Thus the behavior 
described under Phase 1 and Phase 2 is correct only for that portion of 
the deformation that does not generate a midplane force such that 

p 

is greater than one. In terms of Z this condition states that Z 
is less than one, and hence the deflection W must be less than 2 
(see eq. (106)). 

Consequently, equation (117) can be used to calculate the final 
deformation of the plate only if the impulse, I, is small enough so 
that the deformation of the plate is less than* 2. When the plate is 
subjected to impulses greater than this value, a Phase 5 motion will 
begin. The plate will be separated into two parts. In the central 
region where the deflection is greater than two, no bending moments 
are present and the plate behaves as a membrane. In the outer region 
a coupling of midplane forces and bending moments occurs and the plate 
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in this region behaves as it did in Phase 1. The point of separation 
of these regions, | = n, is defined by 

W(n,-r) = 2 (118) 


This point then forms a membrane hinge. This portion of the motion of 
the plate is referred to as Phase 3* 

The time, t 2 , at which the Phase 3 motion begins is defined as 

W(p,x 2 ) = 2 (119) 

Using the expression for Z given in equation (106) and the value 
of the deflection at the membrane hinge, n, given by equation (119 ) > 
the governing equation for Phase 3 ( -q • (89 )) becomes 


6l [ r W__ C d* = -(1 + 3n) - £ [ (120) 

J nL^o J * J n 


In order to evaluate the left side of equation (120) the accelera- 
tion of the outer region (n < | < l) must be evaluated. The velocity 
of the region is given in equation (90) in terms of the position of 
hinge and the transverse velocity at the hinge point W T (n,-r). 

Since the displacement at the hinge is by definition equal to a 
constant 2 (see eq. (ll8)), the transverse velocity is dependent on the 
slope of the deflection curve at the hinge and the velocity of the hinge, 
n T . In addition, as the deformation of the plate is monotonically 
increasing, the membrane hinge must be traveling toward the support. 


Thus 
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W T (n,t) = — n T (121) 

r 1 - n 


With equations (90) and (121) the velocity of the plate in the 


outer region is given by 


W = ■ 2 ~— ■ (1 - |) n < C < 1 (122) 

T (1 - n) 2 


The acceleration can now be written as 


W„ = 


L(1 - n)3 


{(1 - n)n TT + 2n T 2 } 


(1 - 5 ) 


n < 5 < 1 


(123) 


Substitution of equation (123) into equation (120) and performing 
the necessary integration results in the following differential 
equation 


r 2, 

iRl - n)n TT + 2n T 2 ](l + 3n) = -(l + 3n) - l / W^ (124) 

J n 

Equation (124) can now be approximated with the assumption of 
conical surfaces in which W(n,-r) = 2 and W(1 ,t) = 0, to give 


l[(l - n)n TT + 2n T 2 ](l + 3n) = -2(1 - n) 


(125) 
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The solution to this differential equation results in the velocity 
and position of the membrane hinge. The integrating factor for the 
above equation is 


2n x 



(1 

+ nr(l + 3n) 

Multiplying equation (125) by this factor gives 


2 


T ± 

n r 

4 

1 dx 

(1 - n) k 

(1 + 3n)(l - 


dn 
X dT 


Hence the integral of equation (126) is 


(126) 


s? 1 

1 

(1 - n ) 4 

" 192 


A°8_ + + + 81 log 


\{l - n) (i . n ) 2 (l - n)5 


l±3r] 
l-n J 


+ C, 


(127) 


If the Phase 3 behavior begins while the plate is in Phase 2, 
the membrane hinge occurs at the center of the plate, that is, at 
n = 0. The hinge velocity at that time is (from eq. (121)) 


(»T > 0 - \ (128) 

The constant of Integration, Cj, in equation (127) can now be 
evaluated and is given by 

°3 " * I55 (244) 


(129) 
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Substituting equation (129) into equation (127) gives 


n-t 


(1 - n)‘ 


= -F(n) + l(n T ) 2 


( 150 ) 


where 


F(n) = 


1 

108 

192 

(1 - n) 


J2. 


64 


(1 - n) _ n )2 ^ . n )5 


+ 8l log 




244 


(151) 


Values of this function for a range of n from 0 to 1.0 are given 
in appendix C. The motion of the plate stops when n T = 0, therefore, 
the final position of the membrane hinge is given by 

F(np) = l(nr)l (132a) 

or with equation (128) 

F(n F ) = £ A t 2 (t 2 ) (152b) 

The velocity, A^iv,), can be evaluated from equation (ll6), since 
the deflection A is known to be 2 at t = t 2 . Therefore 


A t 2 ( t 2 ) = * + [(St-l) {l2 + (2t x ) 2 } - 32 ]^ (135) 

Equation (152) applies to the condition where the membrane behavior 
is initiated during the Phase 2 motion of the plate. If the impulse is 
great enough, membrane behavior can be initiated during Phase 1, that is. 



64 


before the bending hinge reaches the center of the plate. Since membrane 
behavior occurs when deflections are greater than 2 (see eq. ( 118 ) ) and 
from equation (107) 

W(|,t) = 2t 0 < l < p (134) 

the Phase 3 behavior of the plate is initiated directly from Phase 1 
for all impulses where t becomes greater than unity. 

For this case the midplane forces Eire constant in the central region 
(0 < | < p). The membrane hinge must thus originate at the point 

‘ = P 1 = V 

The differential equation for Phase 3 in this case is identical to 
the previous case. The solution to the differential equation is given 
by equation (127). The initial conditions, however, are different. The 
initial position of the membrane hinge, n e , is equal to the position of 
the bending hinge, p^, at time, t = 1. Thus 

n Q = p(l) = p-l (135) 

The transverse velocity of the plate at the hinge is 2. Thus from 
equation (121) 


(n x ) Q - (1 - r^) - (1 - Pl ) (136) 

Evaluation of the constant of equation (127) for the above condition 


gives 
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( n T )q 
(1 - no) 4 


1 

108 

192 

1 - no 


+ —IZ 

(1 - n Q ) 2 


+ 


64 

(1 - n D ) 5 


+ 8l log 


'l ± 5no\ 

^ - n o / 


(157) 


With equation (137) and noting that n Q = the expression for 
the velocity can be written as 


n T * 

— T 

" <"4 1 

ij. 

— 1 

ij. 

(1 - n) 

-J 


(1 " P 2 ) 


- F(n) + F( Pl ) 


(138) 


where F(n) is defined by equation (ljl). 

The motion stops when n T is zero, hence the position of the 
hinge when the plate comes to rest, np, is given by 


F(np) = I 


( R t)q 

(1 - p -^ 4 


+ F( Pl ) 


(139) 


With the use of equation (136) 


F(np) - 


= I 


(1 - Pl r 


+ F( Pl ) 


(l4o) 


The final position of the membrane hinge is now known. If the center 
portion of the plate does not turn into a membrane until the plate is 
in Phase 2, equations (132) and (133) apply. If the center turns into 
a membrane during the first phase of the motion, equation (l4o) applies. 
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The inner portion of the plate (0 < | < n) during this time behaves 
like a membrane. That is, it carries no bending and transmits only a 
constant tensile load in the direction of its midplane. A very detailed 
and thorough discussion of such behavior is given by Hudson (ref. l) 
and by Frederick (ref. 2). Some of the more pertinent aspects of the 
motion are repeated herein. 

When Phase 3 motion is initiated during Phase 1 motion, the membrane 
portion of the plate behaves as follows. The bending hinge that 
originated during Phase 1 continues to travel inward toward the center 
of the plate. As the hinge sweeps over each annular material element 
in the flat central region, it tilts it into the shape of an annular 
truncated conical element behind the hinge. As the material cannot 
support bending, it is assumed no significant amount of work is done in 
bending. In addition, it is assumed that no thinning of the plate takes 
place. The portion of the plate passed over by the hinge is assumed to 
move in a manner that makes its slope agree with the slope of the outer 
portion of the plate. Thus when the motion of the membrane portion of 
the plate stops, which occurs after the hinge reaches the center, that 
portion forms a cone with a slope equal to the slope of the outer 
portion. 

This conical behavior is also assumed if the Phase 3 motion is 
initiated from Phase 2 motion. Thus 

Ap - 2 2 

np 1 - np 


( 141 ) 
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or 



( 142 ) 



XI. DISCUSSION OF NUMERICAL RESULTS 


Experimental Results 

Before proceeding with a discussion of the numerical results and 
the all-important comparisons with experimental data. It would seem 
appropriate to first discuss the details of the experiments . All of 
the data used herein were obtained from an extensive experimental 
investigation performed by Florence. The details of the experiments and 
the resulting data were discussed In reference 4 from which much of the 
following discussion Is taken. It Is being repeated here for the sake 
of completeness. 

The experiments were very thorough, covering a wide range of 
impulsive loadings and plates of two different materials: aluminum 

and steel. The aluminum plates were 606I-T6 and the steel plates were 
1018 cold rolled steel. 

A sketch of the experimental setup is shown In figure 11. The 
plates were all nominally 1/4- inch thick and 8-1/2 Inches in diameter. 
The specimens were supported on an annular steel plate with an 8-inch 
inside diameter. The impulse was generated by the detonation of an 
8-inch diameter of sheet explosive (DuPont EL-506D) . The sheet 
explosive was cut to fit the opening in the upper support. A similar 
disk of solid neoprene, l/8 of an inch in thickness, was placed between 
the explosive and the specimen. This served to reduce the high peak 
pressure in the shock wave entering the plate and tended to eliminate 
plastic waves, possible changes of material properties, and spallation. 
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The upper support was used to control the rebound. An outside steel 
ring served as a spacer to prevent the edge of the plate from striking 
the upper support during the rebound. 

Dimensions and properties of the plate used in the experiment are 
given in table I. The yield stress was taken to be the average value 
of static tensile tests of both with and acrots grain specimens. B*ch 
experimental stress-strain curve was replaced by a bilinear stress- 
strain curve. The strain-hardening portion of the curves was obtained 
by fitting a straight line through data points above 3- to 4-percent 
strain. 

Twenty-two aluminum and 20 steex plates were tested. The values 
of the initial impulse imparted to the plate were obtained by firing 
free plates in front of a double-flash X-ray unit . The rigid-body 
displacement in the predetermined time between X-ray pictures gave the 
initial plate velocity. 

The permanent central deflections were obtained through the use of 
a traveling microscope. The values of this permanent deflection along 
with the measured initial dimensional impulses are tabulated in table II. 
Also shown in this table are the corresponding values of nondimens ional 
impulse I. 


Large Deflection Effects 

The analyses based on consideration of large deflectional effects 


will be used to predict the behavior of the experimental specimens. 
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TABLE I.- DIMENSIONS AND PROPERTIES OF SPECIMENS* 


Material 

Yield stress, 
lb/in 2 

Denfi i t < y, P 

lb-sec 2 /in 2 

Depth, 

in. 

Radius, 

in. 

Aluminum 

42,000 

0.000253 

0.251 

4.0 

6o6l -T6 





C.R. Steel 

79,000 

0.000732 

0.241 

4.0 

1018 






TABLE II.- EXPERIMENTAL RESULTS* 


Steel plates 

Aluminum plates 

Impulse 

I 

Ebtper . 
Aj> 

nv 0 

lb-sec /m 2 

Impulse 

I 


nV Q 

lb-sec/m 2 

6.97 

8.66 

0.505 

12.71 

13 .42 

0.317 

6.86 

8.43 

0.501 

10.56 

9.94 

O .289 

5.54 

7.44 

0.450 

10.49 

10.96 

0.288 

5.20 

7.14 

0.436 

10.13 

10.61 

0.283 

4.69 

7.00 

0.4l4 

7.53 

8.06 

0.244 

3.52 

6.41 

0.359 

7.28 

8.54 

0.240 

3.33 

5.81 

0.349 

7.28 

8.32 

0.24o 

3.24 

5.54 

0.344 

7.28 

8.4l 

0.240 

3.00 

5-05 

0.331 

6.17 

8.06 

0.221 

2.70 

4.48 

0.314 

6.07 

7-75 

0.219 

2.66 

4.75 

0.312 

4.66 

6.34 

0.192 

2.02 

3.78 

0.272 

4.61 

7.08 

0.191 

1.82 

3.22 

0.258 

4.28 

5.99 

0.184 

1.26 

2.56 

0.215 

2.81 

4.94 

0.149 

0.67 

1.06 

0.157 

2.62 

4.84 

0.144 

0.67 

1.03 

0.156 

2.55 

4.05 

0.142 

0.67 

1.20 

0.156 

2.52 

4.69 

0.l4l 

0.64 

1.49 

0.153 

2.44 

4.27 

0.139 

0.4l 

0.80 

0.123 

2.34 

3.89 

0.136 

0.4o 

0.83 

0.121 

1.91 

3.12 

0.123 




1.76 

3.70 ■ 

0.118 

- 1 



1.47 

3.16 

0.108 


*Data taken from reference 4. 
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Examination of the experimental setup 6hovs that there are no 
mechanisms to provide a midplane force prior to Impulse. Thu6 the 
plates must begin their motion in the pure bending mode. All specimens 
must, therefore, start with a deformat tonal behavior similar to the 
Phase 1 behavior discussed herein. 

The governing equation for the motion of the hinge circle during 
Phase 1 is given by equation (109). A solution to this equation is 
given in appendix B. Numerical results for the position of the hinge 
circle at any time r are given for a range of impulses I from 0.2 
to 12. Typical plots of these values are shown in figure 12. 

The discussion in Chapter X pointed out that Phase 1 motion continues 
until one of two possible conditions occurs. Either the hinge circle 
reaches the center of the plate (p = 0) or the time t reaches a value 
of one. From figure 12 it can be seen that the former condition applies 
to the smaller values of impulse (I < 1.5) while the latter occurs for 
the larger Impulses where I is greater than 1.5. 

In this phase of motion (Phase l) the nondlmensional displacement 
of the central portion of the plate is equal to twice the time t. Thus 
for impulses greater than 1.5 the deflection of the central portion 
reaches a value of 2 (z of one) before the hinge circle reaches the 
center. Once this value of deflection is reached the central portion 
of the plate becomes a membrane. (A nondlmensional displacement of 2 
is equivalent to an actual displacement of 2h, that is, the thickness 
of the plate.) The position of the hinge, p^, when the central portion 
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TABLE III— THEORETICAL RESULTS 


I 

p l 

T 1 

®F 


0.2 


0.192 


0.56 

0.3 

... 

0.280 

... 

O.83 

o.k 

— 

O.366 

... 

1.07 

0.5 

— 

0.447 

... 

1.27 

0.6 

— 

0.520 

... 

1.46 

0.7 

— 

0.595 

... 

1.65 

0.8 

— 

O.651 

... 

1.80 

0.9 

— 

0.710 

... 

1.94 

1.0 

— - 

O.768 

0.044 

2.09 

l.l 

— - 

0.820 

0.103 

2.23 

1.2 

— 

0.868 

0.156 

2.37 

1.3 

— 

0.910 

0.210 

2.50 

1.4 

— 

0.953 

0.240 

2.63 

1.5 

— — 

0.992 

0.274 

2.75 

1.6 

0.044 

... 

0.326 

2.97 

1.7 

0.085 



3-17 

1.8 

0.120 

— 

0.403 

3.35 

1.9 

0.151 

— 

0.454 

3-53 

2.0 

0.179 


0.460 

3-71 

4.0 

0.455 

— 

0.694 

6.54 

6.0 

0.564 

... 

0.773 

8.81 

8.0 

0.627 

— 

O.816 

IO.85 

10.0 


— 

0.842 

12.67 

12.0 

0.699 


0.862 

14.50 

l4.o 

0.722 

“ “ “ 

0.875 

16.00 







Hinge circle 



Figure 12.- Motion of bending hinge for various values of 
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initially becomes a membrane is a function of I as shown in figure 13 
(see also table III). 

For impulses less than 1.5 the hinge circle reaches the center of 
the plate before any portion of the plate can become a membrane 
(t < l). The value of the time, t^, at which this occurs is plotted 
as a function of I in figure 1^ (see also table III). 

Subsequent motion of the plate depends on the magnitude of the 
impulse. For the large impulses (I >1.5) the plate goes directly from 
Phase 1 to Riase 3 behavior. The initial location of the membrane is 
given by in table III. The displacement at the hinge point at this 

time, W(n,T 2 ) is two and the velocity, W T (n,T 2 ) is also two (see 
eqs. (llQ) and (4l)). The final position of the membrane hinge is 
calculated from equation (l4o) and its values are tabulated in table III 
for various values of I. It is also plotted as a function of I in 
figure 15* The final plastic deformations of the center of the plates 
are calculated from equation (14-2) and are tabulated in table III. 

For impulses less than 1.5 the plate goes from the Phase 1 tc the 
Phase 2 motion. The initial center deflection for this phase is given 
by 2 t^ and the velocity by 2. The values of are obtained from 

table III for the appropriate values of I. If the impulse on the plate 
is small enough, the motion of the plate will stop before the central 
part of the plate becomes a membrane. The final deformation is then 
calculated from equation (117). This happens for all Impulses smaller 
than 0.95* 



function of impulse. 





function of impulse. 
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If the impulse is in the intermediate range (0 .95 < I < 1 .5 ) > the 
center of the plate turns into a membrane before the motion of the plate 
ends. For this case the fined position of the membrane hinge is calcu- 
lated from equation (132). These values are shown in table III and 
figure 15. The final deformation as calculated from equation (l42) is 
also tabulated in table III. 

A comparison of the calculated deformations at the center of the 
simply supported plate to those obtained experimentally is shown in 
figure 16 for the small and intermediate impulses and in figure 17 for 
the complete range of impulses from 0 to 14. The circles represent the 
experimental data from the 22 aluminum plates; the squares are from the 
data of the 20 steel plates. The solid curve was obtained from the 
present large deformation theory. The dashed curve is based on Wang's 
plastic bending theory. 

The comparisons between the calculated and experimental deformations 
are very good, especially in the small and intermediate impulse range. 

For impulses less than 4-1/2 the agreement between theory and experiment 
is within experimental scatter. As the Impulse increases above 5, the 
calculated deformations become larger than those measured. This devia- 
tion from experiment can be attributed, at lea6t in part, to the omission 
of strain hardening effects. Inclusion of strain hardening effects 
certainly would cause a reduction in calculated permanent deformation. 
Furthermore, these effects would become increasingly important with the 
additional strains caused by increased Impulse. 
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Figure l6.- Central deflection for small and intermediate impulses and 
comparisons with experimental data. 



Bending theory 
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with experimental data. 
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The bending theory of Wang gives poor agreement even for the low 
values of impulse where no membrane behavior is present . As the value 
of the impulse is increased so that portions of the plate turn into a 
membrane, the difference between the bending and large deformation 
theories becomes substantial. This difference is largely due to the 
neglect of the restraints provided by the midplane forces. This effect 
will be further discussed in the next section. 

Lateral Load Effects 

The influence of the lateral loads on the deformational behavior 
of plastic plates can be seen by examining some of the numerical results 
given in appendix A for the case of plate bending with lateral loads . 

The governing equations for Phase 1 motion for this case were given 
by equation (93) • The equation was solved numerically in appendix A for 
various assigned values of Nq . As the values of Nq are assumed to be 
constant in the analyses, u plastic bending moment is present at all 
times. Thus no membrane behavior is considered regardless of the magni- 
tude of the deformation. Even for the limiting case of N 0 = 1, the 
plate does not behave as a membrane as the flow conditions for the 
interaction case is used, that is, e r = 0. Consequently, equation (93) 
is considered to be the governing equation for the motion of the hinge 
even after the deformation becomes greater than 2 (t >1). 

The final deformation of the plate was calculated from equation (105). 
The results are shown in figure 18, where the central deflection of the 
plate is plotted as a function of the uniform lateral loads, N, for 
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Midplane forces , N 


Figure l8.- Influence of midplane forces on central deflection. 
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various values of the Impulse, I. Examination of these curves shows 
that for small impulses where the deformations are also relatively small, 
the midplane forces do not significantly influence the final deformation. 
However, with increasing values of impulse and correspondingly larger 
deformations, the influence of the midplane forces is greatly increased. 

For example, for an impulse of 10 the calculated central deformation 
of the plate is reduced by almost 60 percent with the introduction of a 
midplane force of only 0.2 of the maximum (N = 0.2) . Once the midplane 
forces reach a value in the range of 0.2 to 0.3, further increases of 
this force are not as effective in reducing the central deformation. In 
the case of an impulse of 10 in the above example, an introduction of a 
midplane force of 0 .2 caused a reduction of deformation from 30 to 12 .5 . 
Further increase of the force to its maximum of 1.0, however, reduces 
the deformation only to a value of 7*0. 


XII. COMPARISON OF LARGE DEFORMATION THEORY TO EXPERIMENTAL 


RESULTS AND OTHER METHODS OF ANALYSES 

This section will be devoted to a critical examination of the 
various methods of analyses that have possible application to deformations 
of plastic plates. Conclusions as to the accuracy of these methods are 
based on the deformational results for impulsively loaded plates . These 
conclusions, however, are believed to be valid for other loading condi- 
tions when applied on a deformational basis. 

Five methods of analysis will be discussed. 

1. The large deformation theory developed herein. 

2. The bending theory of Wang from reference 3* 

3. An analysis based on maximum midplane forces. 

4. The membrane theory of Frederick from reference 2. 

3. The midplane Interaction analysis of Jones from reference 5 • 

The calculated central deflections are shown in figure 19 as a 
function of impulse from all five of these analyses. In addition, the 
experimental results from Florence (ref. 4) are plotted for comparison. 

Large Deformation Theory 

The results from this analysis, replotted from figure IT, are shown 
as a solid line in figure 19. Recall that in this analysis the inter- 
action between plastic bending moments and midplane forces are based 
entirely on the deflect ional behavior of the plate. As long as the 
displacements are less than the thickness of the plate (W < 2) this 
interaction is present. However, when the deflections become greater 
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than the thickness, the central portion of the plate becomes a membrane 
and thereafter supports no bending moment . The motion of the bending 
and membrane hinges, which are the governing factors of the final defor- 
mation, is obtained from solutions of appropriate governing equations. 

From figure 19> it can be seen that the comparisons to theory are 
excellent for deformations of the order of 2 to 3 times the thickness 
(W of k to 6). For larger deformations the theory becomes slightly 
conservative. This, however, should be expected in the light of the 
basic assumptions of the theory. The assumption that the plate material 
is a rigid plastic neglects both elasticity and strain-hardening effects . 
Each of these would tend to reduce the calculated final deformation if 
included. In addition, the assumption that the shape of the outer region 
of the plate is conical tends to underestimate the speed of the hinge 
circles and thus increases the calculated deflection. The inclusion of 
all three of these effects would, therefore, bring the theoretical 
calculations into even better agreement with the experiment. 

Bending Theory 

The results from this analysis are shown in figure 19 as a dash — 
double-dot line. In terms of the parameters used herein, the central 
deflection is based on the equation 

a F = 51 (1^3) 

Figure 19 shows that the comparison between the bending theory and 
experiment is extremely poor even for the lower values of permanent 
deformation. Because of this poor comparison it would seen ill advised 


88 


to use this approach when accurate deflections are desired. As the 
method always gives permanent deformations that are above that obtained 
experimentally, the results of a bending analysis car be used to give 
very conservative estimates of permanent deformations. 

It has been previously pointed out that the behavior of the bending 
hinge circle in Phase 1 is directly related to the final deformation of 
the plate. The larger the speed of the bending hinge the smaller the 
permanent deformation. Thus the conservativeness of the bending theory 
is due primarily to a calculated hinge circle speed tha is too low. 

Maximum Midplane Force Analysis 

The results of this analysis are shown in figure 19 as a dot-dashed 
line. This analysis is identical to that for the large deformation 
except for the interaction relationship of moments and midplane forces. 

The interaction relationship used in the midplane force analysis is not 
dependent on deformations. Instead, it is based on the assumption that 
the deflections will be so large that the entire plate will carry the 
maximum midplane force of x^,(N = l). This, therefore, necessitates 
that the bending moment vanish. The final deformation results for this 
case were discussed in a previous section and correspond to the values of 
Ap for N = 1.0 in figure 18 . 

The comparison of the c entral deformation calculated by the maximum 
midplane force assumption to the experimental results show that this 
method of analysis is unconservative, for the large deformations where the 
method should be most applicable. This uncons ervativeness in calculated 
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deformation reflects a calculated hinge circle velocity that is too 
large. The assumption that only midplane forces are present is, there- 
fore, as poor an assumption as one where only bending moments are 
permitted. 


Membrane Theory of Frederick 

The results from this analysis are shown in figure 19 as a short 
dashed line. In terms of the parameters used herein the central 


deflection as calculated in reference 2 is given by 



The comparison of the experimental data and the theoretical results 
shows excellent agreement for the larger deformations. For small defor- 
mations the agreement is very poor. However, it should be noted that the 
theoretical model used by Frederick was never intended to be used for 
small deformations. 

As the analysis has many assumptions it is difficult for the author 
to fully understand the theoretical reasons for the good agreement 
obtained for larger Impulses as opposed to the very poor agreement in 
the lower impulse range. 

The analysis doe6 assume that only membrane forces are present 
which is the same assumption used by the maximum midplane force 
analysis just described. The behavior of the membrane is also assumed 
to be similar to the Phase 1 and Phase 2 behavior described herein. 
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The outer portion of the membrane in Phase 1 and the total membrane in 
Phase 2, however, are assumed to be rigid in Frederick's approach. Both 
of these assumptions should cause the membrane analysis to underestimate 
the final deformation. The third important difference of the two methods 
is that the membrane analysis assumes that the hinge velocity is uniform 
and dependent only on the initial velocity and the yield stress and 
density of the plate material. Thus it must be this assumption that 
compensates for the unconservatism of the other two. 

Interaction Analysis of Jones 

The results of Jones' analysis are shown in figure 19 as a long 
dashed curve. The data for this curve were taken from a figure in 
reference 5 • 

As can be seen from figure 19, Jones' analysis gives slightly 
unconservative results in the higher impulse or deformation range. 
Addition of the neglected strain hardening an^ elasticity effects would 
result :*r. even more unconservative results. His calculations for small 
permanent deformations are in very pcor agreement . 

Jones assumed in his analysis that since the deformations were 
large, the midplane forces would always be at their maximum. He also 
permitted the plate to carry bending moments during the Phase 1 motion. 
These moments ranged from the maximum M Q at the bending hinge to zero 
at the support. Therefore, he violated the interaction relationship 
based on the Tresca yield condition. His flow conditions were the same 
as those used in the constant midplane-force analysis. Thus it would 
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seem that his calculated deformations should be smaller than those 
showed by the dot-dashed line. The fact that they are not smaller is 
due to an additional assumption used by Jones concerning the hinge 
velocity. In particular, he assumed the hinge velocity to be the 6sme 
as that obtained by Wang for the bending case. As was shown in the 
section dealing with the comparison of the bending theory, this predicted 
hinge velocity is very much smaller than it should be. This lower value 
of velocity must then be responsible for increasing the value of his 
calculated deformation above that of the dot-dashed curve. 



XIII. CONCLUDING REMARKS 


An analysis of a simply supported circular plate under an impulsive 
loading has been presented. The analysis considers the influence of 
both bending and midplane forces and includes large deformation effects. 
Shear and rotary inertia effects are neglected. The plate material is 
assumed to be isotropic and a rigid plastic. The yield criteria ai*e 
based on the Tresca condition. 

Governing equations were developed and solved for three phases of 
motion. The initial Phase 1 included a bending hinge that traveled from 
the support to the center of the plate. Phase 2, which was initiated 
when the hinge reached the center, continued until either the motion 
ceased or a portion of the plate became a membrane. Phase 3 described 
the motion of the membrane hinge from its initial point to its final 
resting position. 

Two types of bending-moment— midplane-force interaction relationships 
were studied: One was based on displacement of the neutral surface from 

the midplane surface. The second was based on the magnitude of midplane 
forces only. The former relationship took account of large deformations 
and midplane distortions and hence allowed regions of the plate to 
behave as a membrane. The latter only considered the lateral forces 
and no membrane behavior was permitted. 

Comparisons of the analytical results with experimental data resulted 
in the following conclusions. 
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1. Interaction relationships based on large deformations should 

be used In order to get reliable data for a complete range of deflections 
and Impulses. 

2 . Analyses of plastic plates based on bending behavior alone give 
highly conservative results and hence should only be used to get con- 
servative estimates of the permanent deformations. 

3. Even a small amount of mldplane force has a significant effect 
on the final deformations of plastic plates. 

4 . Analyses of plastic plates based on maximum membrane forces 
throughout the plate tend to be seriously unconservative. Thi6 method 
should not be used unless accompanied by a conservative assumption to 
counteract this influence. Under no conditions should this assumption 
be used to study deformatlonal behaviors of plates with deform*.. "is of 
the order of one to two plate thicknesses. 

5. For final deformations of the order of five thicknesses and 
greater the influence of strain hardening and elasticity bee canes 
important . 

Elasticity effects can be approximately accounted for in the 
present analyses by incorporating an elastic phase of motion prior to 
the Phase 1 motion considered herein. This elastic phase could be used 
to evaluate the energy dissipated in elastic deformational behavior and 
correct the initial velocity conditions on Phase 1, accordingly. 

A more accurate and sophisticated approach, however, would be to 
permit the plastic behavior to initiate at the center of the plate and 



then spread toward the supports. To do this it would be necessary to 
define another hinge circle separating the elastic and plastic portions 
of the plate. 

Another possible area of refinement is in the analysis of the 
membrane portion of the plate. An initial effort in this area should 


deal with the influence of reduced thickness . 
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XVI . APPENDIX A 



The governing equation for Phase 1 in terms of the impulse I and 
the lateral load, N 0 , is (see eq. (93) 

_ -(1 - N e 2 ) - 2 N q t( 1 + p) 1 

Pt= I(l-p)(l+3p) 


The initial condition for this phase is that p = 1 for r = 0. 
This phase ends at t = when p = 0. Phase 2 motion then begins 

and continues until a final deformation is reached. The equation for 
the final deformation is given by equation (105) as 

N Q Ap 2 + 2(1 - N Q 2 )Ap = 21 + 4^(1 - N 0 2 ) + 4 t 1 2 N 0 (A-2) 

For the case of pure bending (N 0 = 0) the solutions for 
equation (A-l) can be found in closed form as 


l_p_p2 + p 3=i 


(A-3) 


so that 



(A-4) 
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100 


The solution of equation (A-2) for Ng = 0 is then 

Ap = 31 (A-5) 

For the case of the maximum midplane force (Ng =1) the solution of 
equation (A-l) is 

t 2 = l|^(l + p) 2 - 8(1 + p) + ^ log - — g — + 10] (A-6) 

so that 

t x 2 = 0.727^1 (A-7) 

The solution of equation (A-2) then can be written as 

Ap 2 = 21 + ^r^ 2 (A-8) 

A p = 2.21 \/f (A-9) 

For the general case equation (A-l) is solved numerically for 
and the value of A p calculated from equation (A-2) . The flow 
diagrams, computer program, and numerical results follow. 
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000075 

000076 

000077 

000102 

000106 

000115 

000115 

000117 

000123 

000125 

000130 

000130 

000133 

000156 

000172 

000172 

000174 

000176 

000200 

000202 


COEFS ( 2 ) *- 1 • 0 

C0EFS<3>»-1.0 

COEFS ( 4) »1 • O-OTAU/F 1 

CALL FALG ( COE FS 1 3 , 0* R , TEMP , I ERR ) 

IF ( IERR.NE *01 PRINT 5« I ERR 

5 FORMAT ( //I OX*ERROR IN FALG SUBROUTI NE*/ ) 

00 6 M-1,3 

IF (AIMAGIRIM) I.EO.O.O) GO TO 8 

6 CONTINUE 
PRINT 7 

7 FORMAT I //10X*ALL ROOTS COMPLEX*/) 

8 RHO(l) -REAL (RIM)) 

CALCULATE INITIAL RHOT VALUE 

RHOT(l)— l2.*NTHETA*TAUCll*(l.*RHOU))m.-N1HETA**2))/IFI*(l.- 
1RH0U) )*C1 • *3. *RHO( 1 ) ) ) 

PRINT 9, Nt TAU(N) »RHQ(N) tRHOT (N) 

9 FORMAT (4X*N*» 1 IX*TAU*,17X*RH0*, I7X*RHOT*//2XI3 t 3! 5XE15. 8) ) 

C BEGIN ITERATION ON RHO. TAU» RHOT 

10 N-N+l 

0ELRH0«M0T(N-1)*0TAU 
RHO (N)-RHO! N— 1 ) ♦DEL RHO 
TAU ( N) -TAU ( N-I >*DTAU 

RHOT INI— ( 2.*NTHETA*TAU(N)*( 1«*RH0(N) )♦( l.-NTHETA**2) ) / 
l (FI*( 1 •— RHO(N) )*! !• ♦3.*RH0C N) ) ) 


C OUTPUT 

* RHO -RADIUS OF HINGE CIRCLE 

* RHOT - VELOCITY OF HINGE CIRCLE 

* TAU - TIME OF HINGE CIRCLE 


000227 

000247 


000247 

000253 

000254 

000254 

000260 

000260 

000261 


IF (MOOIN. IPFN1.EQ.0) PRINT li t N.TAU(N) *RHO(N) *RHOT(N) 

11 FORMAT ( 2X13 «3( 5X615.8) ) 

C IF MAXIMUM NUMBER OF ITERATIONS IS REACHED OR RAOIUS RHO) 

C BECOMES NEGATIVE — STOP 

IF (N.GT.MAXN) GO TO 12 
IF (RHO(N).LT.O.O) GO TO 20 
GO TO 10 

12 PRINT 13 

13 FORMAT (//10X*MAXIMUN N REACHED*//) 

GO TO 14 

C INTERPOLATE FOR TIME AT WHICH RHO-O 

20 T AUA-T AU (N- 1 ) *RHO I N— l ) *DT AU/ ( RHO ( N- 1 )— RHO ( N ) I 
CALCULATE AF 


000270 C0(1)=NTHETA 

000272 C0( 2 ) = 2.*C l .-NTHETA**2» 
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-UNSATISFIED EXTERNALS REFERENCES 


Ill 


INFLUENCE OF UNIFORM MID-PLANE FORCES ON 
PLASTIC BENDING UNDER IMPULSIVE LOADING 

ROBINSON— KRUSZENSKI , SKD-A1998, FEBRUARY 1968, RDV-122 
DATE 03/07/68 TIME 19.16.59. 


DT AU= .010 1= 4.0 NTHETA= .10 


N 

TAU 

RHO 

RHOT 
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l • 000000 OOF— 02 

9. 64 3 2505 4E— 01 

-1.789164446+00 

10 

1.000000006-01 

8. 806086 19E— 01 

—5. 90849807E-01 

20 

2* 000 000 OOE— 01 

8.2927601 IE— 01 

—4. 463684 95 E- 01 

30 

3 • 000000 OOfc—O 1 

7. 877619366-01 

—3. 84295030E— 01 

40 

4* QOOOOOOOE— 01 

7. 51 066650E— 01 

— ? . 48865004E— 01 

50 

5 • 00 0000 OOE- 01 

7. 172 861 OOE— 01 

— 3.2 59 34 887 E- 01 

60 

6. 000000 00E-01 

6. 854526 08E— 01 

—3. 100406406-01 

70 

7, 000000006—01 

6. 549939 72E— 01 

— 2.9857645 IE- 01 

80 

8. OOOOOOOOE-Ol 

6. 25 53 731 4E— 01 

— 2.9012 828 OE- 01 

90 

9. OOOOOOOOE-Ol 

5# 96 82205 7E— 01 

—2 . 83 85 9153 E— 01 

100 

i.OOOOOOOOE+OO 

5.68655 887E- 01 

-2* 792425 42 E— 01 

110 

1.1000Q000E+G0 

5.408902 56E— 01 

—2. 75932 814 E— 01 

120 

1 • 200000 OOE+OO 

5. 134057 056—01 

—2*73696801 E- 01 
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l. 30000000E+00 

4* 86 102 52 7E— 01 

—2. 723753246—01 
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4. 588944 7 5E- 01 

-2.71860561E-01 

150 
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—2. 720823926-' 1 
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1 .600000006+00 

4. 044604 08E— 01 

-2.730001546 )1 
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1.70000000E+00 

3.77 094097E— 01 

-2.74597925E 01 

180 

1.80000000E+00 

3. 495372 5 5E— 01 

-2. 76 882340 £-01 

1 v 0 

1.9 00000 OOF +00 

3.21 720101E— OL 

-2. 798824806-01 

200 

2 • OQOOOOOOE+OO 

2.93 568 86 9E— 01 

-2.836518216-01 

210 

2 . ldooooooE+oo 

2* 650032 09E-01 

-2.882723806-01 

220 

2.20000000E+00 

2.35933030E— 01 

-2.938620616-01 

230 

2. 30000000E+00 

2. 062 5442 8E- 01 

-3.00566357E-01 

240 

2.4 3 00 00 OOE +00 

1 • 75844160E— 01 

-3.086771316-01 

250 

2.50000000E+00 

1 • 4455 1752E- 01 

-3.184632066-01 

260 

2 . 600000 OOE + 00 

1.1218771 IE— 01 

-3.304219856-01 

270 

2. 70000000E+00 

7. 85050049E— 02 

-3.452709346-01 

280 

2.80000000E+00 

4. 31 632146E- 02 

— 3. 641406686—01 

290 

2 . 900000 OOE+OO 

5 • 6984 1 823E— 03 

— 3.889 3 156 9E-01 


T AUA= 2.9 1 46 1 683 E ♦00 AF ( l ) = 8. 1944409IE+00 0. 
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.010 I- 6.0 

NTHET A= .10 


H 

TAU 

RHO 

RHOT 

1 

1. 00000000E-02 
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10 
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?0 

2. 000000 OOE-Ci 
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-3. 57967732 E- 01 

30 

3. OOOOOOOOE-OI 

8 • 28546299E— 01 

-3. 06 689 398 E- 0 1 

40 

4.00 0 000 00 1 - 0 1 

/ . 99 320530E— 0 1 

-2.77152995E-01 

50 

5. OOOOOOOOE-Oi 

7* 72 536510E— 01 

-2. 57796503 E- 01 

60 

6. OOOOQQOQF— 01 

7.47407232E— 01 

—2.44 148 3 3 OE- 01 

70 

7. OOOOOOOOE-OI 

7* 23 469 545 E- 01 

-2. 340 71 95 8E— 01 

80 

3.00000000F-01 

7 . 004239 16E— 01 

-2. 26404202E— 0 1 

90 

9.00000000t-01 

6.78063379E-01 

-2. 204526 73 E- 01 

100 

l.OOOOOOOOE+OO 

6.5623764 IE-01 

-2. 15 7 77627 E- 01 

110 

1. 10009000E+00 
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-2. 1208578 7 E- 01 

120 

1*2 0 00000 OE +00 
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1* 3000 00 OOF + 00 
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1 • 50000000E + 00 
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160 

l.oOOOOOOOE+OO 
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4.09973721E-01 

—2. 0399453 3E—0 1 
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2. 300000 COE+OO 
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— 2.051 04 988 E- 01 
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-2.0811 8830E-01 
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2 • 600000 OOE + 00 

3.27 3409 63E— 0 1 
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270 
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290 
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-2. 46 1710 03 E- 01 
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3. 60000000E+00 

1. 00 4473 34E— 01 

-2. 5 3 7563 3 6 E- 01 
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3. 70000000E+00 
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-2. 62 701583 E- 01 
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3. 8 00000 OOE +00 
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-2* 7 34076 01 E- 01 
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3. 9000000 OE +00 

2.00405872E-02 

-2. 86 46 844 3E— 01 


T AUA= 3. 96887600E *00 AF<1I= 1 . 10328784E + 01 0. 
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DTAU= 

• 010 1= 8.0 

NTHETA* .10 
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RHO 

RHOT 
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1.00000000E-02 

9.748412 59F- 01 

-l.?5834'55E+00 

10 

1.00000000E-01 

9.16326172*E-0l 

-4.U9768396E-01 

20 

2. OOOOOOOOE-Ol 

8. 80865910E— 01 

-3. 06 837213 E- 01 

30 

3. 000000 00E-01 

8.524241 19E-01 

—2. 62 19 366 7E- 01 

40 

4. OOOOOOOOE-Ol 

8.27465145E-01 

-2. 363 78 952 E- 01 

50 

5. OOOOOOOOE-Ol 

8. 04 6445 2 7£- 01 

-2. 19375231E-01 

60 

6* OOOOOOOOE— 01 

7.83281360E-01 

-2. 0730672 5E- 01 

70 

7. OOOOOOOOE-Ol 

7.62 975496E-01 

-1. 98 32063 IE- 01 

80 

8. OOOOOOOOE-Ol 

7. 434687 16E— 01 

— 1.9 140742 OE- 01 

90 

9. 000000 00E-01 

7.245831 80E- 01 

—1. 8596542 3E- 01 

100 

1.00000000E+00 

7. 06 190328E-01 

-1. 81611781 E- 01 

110 

1. 10000000E+00 

6. 88 193 65 OE— 01 

—1. 7809041 IE- 01 

120 

l « 20 0000 OOE+OO 

6.70518442E-01 

- 1. 752231 93E- 01 

130 

1.30000000E+00 

6.531053 76E-01 

-1.72882294E-01 

140 

1.40000000E+00 

6. 35 90628 IE- 01 

—1.709 736 17E- 01 

150 

1 • 500000 OOE+OO 

6.18881280E-01 

— 1. 6942 6458 E- 01 

160 

1 • 60 00000 0E+ 00 

6. 01 996768E— 01 

-1.68186795E-01 

170 

1. 70000000E+00 

5. 85 22 396 2E— 01 

-1.672128 15 E- 01 

180 

1 • 600000 OOE+OO 

5.68537832E-01 

— 1. 66471834E-01 

190 

1 • 900000 OOE+OO 

5. 51 9162 84E— 01 

- 1.659381 40E- 01 

200 

2. 00000G00E+00 

5. 35 339542E-01 

— 1*65591453 E- 01 

210 

2 • ICO 000 OOE+OO 

5.18789656E-01 

-1.6541581 IE- 01 

220 

2.20000000E+00 

5.022501 07E-01 

— 1.65398746E-01 

230 

2.30000000E+00 

4. 85 705492E-01 

- 1.655306 85E- 01 

240 

2. 400000 OOE+OO 

4.69141248E-01 

— 1.65 8044 99E- 01 

250 

2. 50000000E+00 

4. 52 54343 IE— 01 

- 1. 662 1 51 80£- 01 

260 

2.60000000E+00 

4. 35 8985 09E— 01 

-1. 66759608E-01 

270 

2.70000000E+00 

4. 19193184E— 01 

-1.67436390E-01 

280 

2. 800000 OOE+OO 

4.02414222E-01 

-1.68245765E-01 

290 

2. 900 000 OOE+OO 

3. 85 5482 95 E- 01 

— 1 .691 89566E-01 

300 

3 • OOOOOOOOE *-00 

3. 68581 8 2 OE- 01 

- 1. 70271225 E- 01 

310 

3. 1 OOOOOOOE+OO 

3.51 500794E-01 

— 1. 7149584 IE- 01 

320 

3.20000000E+00 

3.342906?^ '-01 

— 1. 72870294E-01 

330 

3. 3 00000 OOE+OO 

3. 16935 9o0E— 01 

- 1. 7440342 IE- 01 

340 

3 .400000 OOE+OO 

2.99420442E 01 

— 1.76 1 06262 E- 01 

350 

3. 50 0000 OOE+OO 

2.61726516E-01 

-1. 779923 94 E- 01 

360 

3 • 600000 OOE+OO 

2.63835133E-01 

-1.80078373E-01 

370 

3* 70000000E + 00 

2.45 725437E-01 

- 1. 82 38431 OE- 01 

380 

3. 80000000E+00 

2. 27 374371 E- 01 

— 1.84934637E-01 

390 

3. 900000 OOE+OO 

2.08 756207E-01 

—1.877 59 11 8E-01 

400 

4 • OOOOOOOOE +00 

1 • 89841949E— 0 1 

— 1 • 908 942 02 E- 01 

410 

4. 10000000E+00 

1.70598595E-01 

-1.94 3 848 5 IE- 01 

420 

4. 2 00000 OOE+OO 

1 .50988 17 2E- 01 

- 1.9828706 IE- 01 

430 

4. 30000000E+0G 

1. 30966487E-01 

-2.02671374E-01 

440 

4.40000000E+00 

1. 1048146 IE-01 

-2.07627906E-01 

450 

4.50 0000 OOE+OO 

8.94 708677E-02 

-2.13273698E-01 

460 

4. 60 0000 OOE+OO 

6.78591 765E-02 

— 2. 19 763777E-01 

470 

4.70000000E+00 

4. 555304 86E-02 

— 2# 273083 5 5E- 01 

480 

4. 800C0000E+00 

2.24346799E-02 

-2.36200640E-01 

490 

4. 900000 OOE+OO 

-1.64840441E-03 

— 2.46863946E-01 


T AUA* 893291 04E *00 AF (II* l . 35000324E+01 0. 
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OTAU* .010 I* 10.0 NTHETA- .10 


N 

TAU 

RHO 

RHOT 

l 

l • 000000 OOE -02 

9.775125416-01 

-1.123966776*00 

10 

l.OOOOOOOOE-Oi 

9. 25325755E-01 

-3. 64 759686E-01 

20 

2.00000000E-01 

8.937896046-01 

-2.72570989E-01 

30 

3.00000000E-01 

8.685431626-01 

-2.32520735E-01 

40 

4.000000006-01 

8.46423763E-01 

-2.093125126-01 

50 

5.00000000E-01 

8.262288956-01 

-1.939838126-01 

60 

6.000000006-01 

8. 07 34 973 6E- 01 

-1.830665266-01 

70 

7. 00000000E-01 

7.894 o712E-0l 

-1.749023516-01 

80 

8.000000006-01 

7.722352966-01 

-1.685874436-01 

90 

9.00000000E-01 

7.556109356-01 

-1.635828666-01 

100 

i.ooooooooe*oo 

7.39441233E-01 

-1.595454506-01 

110 

1. 10000000E*00 

7.236405296-01 

-1.562454316-01 

120 

1 • 20000000E *00 

7.081427406-01 

-1.535227686-01 

130 

1.300 0000 OE *00 

6.92 895607E-01 

-1.512623756-01 

140 

1 • 400000006*00 

6.778569306-01 

-1.493792986-01 

150 

1.5 00 0000 06 *00 

6.629920056-01 

-1.478094486-01 

160 

1 • 600000006*00 

6.48271829E-01 

-1.465035966-01 

170 

1. 700000006*00 

6.336718076-01 

-1.454233466-01 

180 

1.800000006*00 

6.191 70799E— 01 

-1.445383726-01 

190 

1.900000006 4*00 

6. 04 7504 1 2 E— 01 

-1.438244766-01 

200 

2 • 000000 OOE *00 

3. 9039444 8E— 01 

-1.432621906-01 

210 

2. 100 OOOOOE 4*00 

5. 76 08 8482 E— 01 

-1.428357576-01 

220 

2 • 20000000E 4-00 

5. 61 819528E-01 

-1.42532377 E- 01 

230 

2.300000006*00 

5.475757656-01 

-1.423416396-01 

240 

2.400000006*00 

5.333463256-01 

-1.422550896-01 

250 

2.500000006 4*00 

5. 19 1 21 1 08E-0 1 

-1.422659056-01 

260 

2 . 6 00 000 OOE 4* 00 

5.048906346-01 

-1.423686416-01 

270 

2 • 700000 OOE 4*00 

4.906459106-01 

-1.425590326-01 

280 

2. 800000006*00 

4. 76378326C-01 

-1.428338476-01 

290 

2.900000006 4*00 

4.620795506-01 

-1.431907716-01 

300 

3.000000006*00 

4.477414446-01 

-1.436283256-01 

310 

3. 10000000E*00 

4. 3335 59 8 3E -01 

-1.441458036-01 

320 

3.200000006 4*00 

4.18915180E-01 

-1.447432356-01 

330 

3.300000006*00 

4.044110126-01 

-1.454213626-01 

340 

3.400000006*00 

3. 89 83 534 9E- 01 

-1.461816366-01 

350 

3. 5000000oE*00 

3.751798856-01 

-1.470262356-01 

360 

3 • 60 0000 OOE *00 

3.60436062E-01 

-1.479580866-01 

370 

3. 700000006*00 

3.45594994E-01 

-1.489839206-01 

380 

3.800000006*00 

3.30647383E-01 

-1.500993376-01 

390 

3.900000006*00 

3.155834296-01 

-1.513188986-01 

400 

4.000000006*00 

3.00392729E-01 

-1.526462406-01 

410 

4. 10000000E*0C 

2.850641606-01 

-1.540892296-01 

420 

4.200000006*00 

2. 69 585754 E- 01 

-1.556571456-01 

430 

4.300000006*00 

2.539445406-01 

-1.573609286-01 

440 

4.400000006*00 

2. 38 1263 70E-01 

-1.592134716-01 

450 

4.500000006*00 

2.221157016-01 

-1.612300136-01 

460 

4.600000006*00 

2. 0589534 7E-01 

-1.634286236-01 

470 

4. 700000006*00 

1.894461646-01 

-1.65830838E-01 

480 

4.800000006*00 

1.727466696-01 

-1.684624886-01 

490 

4.900000006*00 

1.55772^626-01 

-1.713547826-01 

500 

5. OOOOOOOOE^OO 

1.384961326-01 

-1.745457676-01 

510 

5. 100000006*00 

l. 20 88 548 3E- 01 

-1.780823006-01 

520 

5.200000006*00 

1.029035426-01 

- 1 • 82022784E-01 

530 

5.300000006*00 

8.450673816-02 

-1.864410146-01 

540 

5.400000006*00 

6.564323676-02 

-1.914317316-01 

550 

5.500000006*00 

4.625049686-02 

-1.971188426-01 

560 

5.600000006*00 

2. 62 51 83 80E-02 

-2.036679686-01 

570 

5. 700000006*00 

5.551453056-03 

-2.113062 8 BE- 01 


T AUA= 5. 7 26 1 8326E *00 *F I 1 > - l . 571 1 00A6E+01 0. 
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DT X\J~ .010 I* 12.0 NTHETA- .10 

N TAU RHO RHOT 

l L.00000000E-02 9.79482068E-01 - 1. 02501 1416*00 


10 

l.OOOOOOOOE-Ol 

9. 

20 

2. 000000 OOE -01 

9. 

30 

3. 000000 OOE-Ol 

8. 

40 

4. OOOOQOOOE-Ol 

8. 

50 

5. 00 0000 OOE -01 

8. 

60 

6. OOOOOOOOE-Ol 

8. 

70 

7. OOOOOOOOE-Ol 

8. 

80 

8. OOOOOOOOE-Ol 

7. 

90 

9. OOOOOOOOE-Ol 

7. 

100 

l. 000000006*00 

7. 

no 

1. 100000006*00 

7. 

120 

1 • 2 000000 OE *00 

7. 

130 

1. 300000006*00 

7. 

140 

1. 400000006*00 

7. 

150 

1 . 50000000E 4-00 

6. 

160 

1.600000006*00 

6. 

170 

1 • 70000000E *00 

6. 

180 

l .80000000E*00 

6. 

190 

1 • 90000000E*00 

6. 

200 

2 • 000 000 OOE +00 

6. 

210 

2.10000000E+00 

6. 

220 

2 . 20 00000 OE* 00 

6. 

230 

2.300000006*00 

5. 

240 

2 .400 000 OOE 4-00 

5. 

250 

2.50000000E4-00 

5. 

260 

2 * 600000 OOE *00 

5. 

270 

2. 70000000E*00 

5. 

280 

2 • 800000 OOE 4- 00 

5. 

290 

2. 900000006*00 

5. 

300 

3 • 00 0000 OOE *00 

5. 

310 

3. 100000006*00 

4. 

320 

3 • 2000 J000E*00 

4. 

330 

3.30000000E4-00 

4. 

340 

3.400000006*00 

4. 

350 

3 • 50 00 00 OOE 4-00 

4. 

360 

3. 6000000064-00 

4. 

370 

3. 70000000E4 , 00 

4. 

380 

3.80000000E*00 

4. 

390 

3* 900000 006*00 

3. 

400 

4. 00 0000 OOE ♦ 00 

3. 

410 

4. 10000000E4-00 

3. 

420 

4. 2 00000 OOE tOO 

3. 

430 

4. 30 0000 OOE 4-00 

3. 

440 

4. 4 00000 OOE 4-00 

3. 

450 

4. 500000006*00 

3. 

460 

4. 60 0000 OOE ♦ 00 

2. 

470 

4. 700000006*00 

2. 

480 

4. 80000000E*00 

2. 

490 

4. 90 0000 OOE 4-00 

2. 

500 

5 • 00 0000 OOE 4- 00 

2. 

510 

,100000006*00 

2. 

520 

5 • 2 0000000E *00 

2. 

530 

5.3 00000 OOE 4* 00 

1 • 

540 

5. •* 00000 OOE *00 

X * 

550 

5.50000000E+00 

L • 

560 

5.60 0000 OOE + 00 

1 • 

570 

5 • 7 00000 OOE *00 

1 • 

580 

5 • 80 0000 OOE 4- 00 

X • 

590 

5.900000006*00 

1 • 

600 

6. 000000006*00 

8. 

610 

6. 1 OOOQOOOE+OO 

6. 

620 

6.20000000E*00 

5. 

630 

6* 3 00000 OOE *00 

3. 

640 

6. 4 00000 OOE *00 

1. 

650 

6. 50 C 000 OOE *00 

-1. 


31942 7596-01 - 3. 318 i 86436-01 
03273630E-01 -2.47589073E-01 
80353384E-01 -2. 10957800E-01 
60294723E-01 -1. 89702092E-01 
419998 75E-01 -1 . 756390 72E- 01 
249131 35E-01 -1 . 6560 19 72E-01 
08708225E-01 -1. 5807638 IE-01 
93174896E-01 - 1 .52236742E-01 
7816865 OE-OL -1.47590687E-01 
63585329E-01 -1 .43824595E-01 
49 347046E-01 -1. 407283 93E-01 
35393813E-01 -l. 381 55668E-01 
21678299E-01 - 1. 36001040E-01 
08162383E-01 -1.341 Q6629E-01 
948148 24E— 01 -1. 3265 3591 E- 01 
816096 17E-01 - l . 3. 35662 7E-01 
68524819E-01 -1.302603068-01 
5554168 IE-01 -1 .293365 38E-01 
42 644003E— 01 - 1 . 28562795E- 01 
29817633E-01 - 1 .27920829E-01 
17050089E-01 - 1. 27395742E-01 
04330254E-01 -l .2697528 7E-01 
916481 35E-01 - 1.2664934 7E-01 
78994672E-01 - 1 .2640953 0E-01 
66361576E— 01 - 1.26248869E-01 
r 3741202E-0l -l .26161575E-01 
4l 1 26435E— 01 - 1.2614285 IE-01 
28510604E-01 - 1 . 26 1 88745E- 01 
15 8873 96E- 01 - 1 . 2629603 3E-01 
03250795E— 01 - 1 .264621 2 IE-01 
90595015E-01 -1.26684978E-01 
77914454E-01 - i . 2696 3074E- 01 
65203638E-01 - l . 27295335E-01 
524571 85E-01 - l . 27681 1 1 OE- 01 
39669756E-01 - 1.28120150E-01 
26836024E-01 - 1 , 28612 584E-01 
1395062 8E-01 - 1.2915891 BE -01 
01008143E-01 -1.29760029E-01 
88003037E-01 - 1. 304171 71E- 01 
74929637E-01 - 1 . 31 1 3 1 988E-01 
61 782086E— 01 -l. 31 90653 3E- 01 
48554303E-01 - l . 32 743297E- 01 
35 2 39937E- 0 1 - 1 . 33645240E- 01 
21 832318E-01 -1.346158416-01 
08324403E-01 - l . 356591 55E- 01 
94 7087156-01 - 1 . 36779378E-01 
809772 75E-01 -l. ^98 3442E- 01 
57 121524E-01 -1.392761126-01 
53 1 32235E-0 l - 1 . 40665 12 36- 01 
38999407E-01 - 1 . 421 5883 90- 01 
247121466-01 - 1.43 76694 7E-01 
102585 19E- 01 - 1.4550071 IE- 01 
956253 88E-01 -1.473732796-01 
807901976-01 - 1. 494000 736- 01 
657607316-01 - l. 515992906-01 
50494805E-01 -1.539925456-01 
349798926-01 -1.566057126-01 
191926446-01 -1. 594700316-01 
03106296E-0I - l .626235966-01 
6689891 16—02 - 1.661 13382E- 01 
990727966-02 - 1. 699980596-01 
27157759S-02 - 1 . 74351 9786- 01 
506435 13E-02 -1.792709636-01 
68911267E-02 -1 . 84880937E-01 
88018629E-03 -1.913512286-01 


TAUA* 6.490138516 *00 


AF ( 1 ) * 1 . 7731631 8E*01 0 
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JTAU = 

•010 1= 4*0 

NTHETA= .30 


N 

TAU 

RHO 

RHOT 

1 

1. OOOOOOGOE-02 

9.643250546-01 

-1.65930084E+00 

10 

1 *000000008-01 

8. 830794556-01 

-5.993991376-01 

20 

2.00000000E-01 

8.298506926-01 

-4.756176186-01 

30 

3.000000008-01 

7.847966846-01 

-4.264110596-01 

40 

4. OOOOOOOOE-Ol 

7.434337016-01 

-4.007136496-01 

50 

5 • OOOOOOOOE-O 1 

7. 0409636 IE— 01 

-3.858102896-01 

60 

6.00000000E-01 

6.659543026-01 

-3.768994826-01 

70 

7 • 000000 00E-01 

6. 285 2 106 IE- 01 

-3.717558216-01 

80 

8,000000008-01 

5.914779816-01 

-3.692173556-01 

90 

9. OOOOOOOOE-O 1 

5.545907576-01 

-3. 686346426-01 

100 

UOOOOOOOOE + OO 

5.177002546-01 

-3.696360886-01 

110 

1. IOOOOO0OE+OO 

4.806403466-01 

-3.72017456E-01 

120 

1.2 GOOOOOOE+OO 

4. 43283950E— 01 

-3.756872376-01 

130 

1 . 3<)ooooooe+oo 

4.05503099E-01 

-3. 806405816-01 

140 

1 .4 GO 000 DOE + 00 

3.671667966-01 

- 3. 869503936- 01 

150 

1,500000 >0E + 0 0 

3.281331546-01 

-•>.947714896-01 

160 

1.60000000E+00 

2.88240532E-01 

-4.043582606-01 

170 

1. 70000000E+00 

2. 472961476-01 

-4. 161007096-01 

180 

1. 80000000E+00 

2. 050598496-01 

-4. 305908226-01 

190 

l • 900000 OOF + 00 

1.612190316-01 

-4.487455836-01 

200 

2. 000000 OOE +00 

1.1 5 3466 93 E— 01 

-4. 720473736-01 

210 

2* 10009G00E + 00 

6. 682523936-02 

-5.030556746-01 

220 

? • 200000 OOE +00 

1.469280626-02 

-5.466393186-01 


T AUA= 2.22664o566*00 AF ( 1 ) •= 6. 061 5 1 369E ♦00 0. 
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0TAU= 

.010 l = 6.0 

NTHE T A= .30 


N 

fAU 

RHO 

RHOT 

l 

i.OOOQOOOOE-02 

9. 709203 U5E-U1 

-1. 35028152E+00 

LO 

1 .OOOOOOOOE-Ol 

°.05042592F-01 

-4.63920921E-01 

20 

?. 000000 001 -01 

8.621621 69E- 3 1 

-3. 821 352566-01 

30 

3. OOOOOOOOE-Ol 

8.2603 02 95E-01 

-3.4118 9 167E-01 

40 

4. 00000000£-01 

7.92992244E-01 

-3. 19363167E-01 

50 

5. OOOOOOOOE-Ol 

7. 616 96653 E- Cl 

-3.06255518E-01 

60 

6. OOOOOOOOE-Ol 

7.31 4 759 09E-0 1 

-2.97926021E-01 

70 

7. OOOOOOOOE-Ol 

7 . 01 9452 58E-0 1 

-2.9253557^E-01 

80 

6. OOOOOOOOE-Ol 

6. 728593 3 3r— 01 

-2.89109391E-01 

90 

9. OOOOOOOOE-Ol 

6. 74 04 938 2 E- 01 

-2. 8708549 IF- 01 

10 0 

l.OOOOOOOOE+QO 

6. 15391936E-01 

-2. 861210; V_ -01 

110 

1. 10000000E+00 

5.86 79 15 96E-0 1 

-2.859993t?5E-0l 

120 

1 • 2000000 0E + 00 

5.58170838E-01 

—2. 865316 88E-0 1 

130 

1 • 3000000 OF + 00 

5.29463555E-01 

-2.87780342E-01 

140 

1.40000000E+00 

5. 006106 99E-01 

— 2.89 54 34 50 E-01 

150 

1.50000' 'OE + OO 

4.71557143E-01 

-2.91846087E-01 

160 

1.60000- jOE+OO 

4 • 42 2492 35E“ 0 1 

-2. 946 8536 OF -01 

170 

1.70000000E+00 

4. 12632710E-01 

-2. 980781 42 E-01 

180 

1 . dOOOOQOCii + OO 

3 • 8265 0744 E-0 1 

-3.02060 86 8E -01 

190 

1 .90000000E+00 

3.52241934E-01 

-3. 06 69 1195 E-01 

200 

2.00000000E+00 

3.21 338031 E-01 

- 3. 12051730 E-01 

210 

2. 10000000E+00 

2.89 861 1 8 IE- 0 1 

-3. 18256398 E-01 

220 

2.2Q000000F+00 

2.57720374E-01 

-3*2 54605 77E- 01 

230 

2.30000000E+00 

2.24806645E-01 

-3.33877083E-01 

240 

?.40000000Ef 00 

1 • 909862 7 IE- 01 

-3.43301 612 E-01 

250 

2 • 500 000 OOF +00 

1.56090749E-01 

-3.55656 06 7E-01 

260 

2.b0000000E+00 

1. 19901286E-01 

-3.7 005 9969 E- 01 

270 

2 • 700090 OOE *00 

8.21235722E-02 

-3.87967205E-01 

280 

2.80000000E+00 

4. 23439742 E- 02 

-4* 10932 380E- 01 

290 

2* 90000000F + 00 

-5* 2938593 5E-05 

-4.416980846-01 


T AUA= 2. 899879186 «-00 AF ( 1 ) = 7. 8305291 8E+00 0. 
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0T AU= 

• 010 1= 8.0 

NTHET A= .30 


N 

TAU 

RHO 

PHOT 

I 

i.ooooooooe-02 

9.74841259E-01 

— 1. 16 7062 06 E> 00 

10 

i.OOOOOOOOE— 01 

9 . 1 8020904E- 01 

-4. 16354C78E- 01 

20 

2.0000G000E-01 

8. p l 1 72035E— 01 

- 3. 279 063526- 01 

30 

3. OOOOOOOOE-Ol 

8. 50 1982 3 OE— 01 

-2. 921296 OOE-Ol 

40 

4.00000000E-01 

8.21936426E-01 

-2. 7288722 9E -01 

50 

5. OOOOOOOOE-Ol 

7.95218934E-01 

- 2. 61 165355 F- 01 

60 

6. OOOOOOOOE-Ol 

7.6947095 IE- 01 

-2. 535465i 66-01 

70 

7. OOOOOOOOE-Ol 

7. 44 3 62 84 9E— 01 

— 2. 4 3 '*327406-01 

80 

8. OOOOOOOOE-Ol 

7. 1968646 IE-01 

— 2. 4*9740776- 01 

90 

9. OOOOOOOOE-Ol 

6.95 30053 IE- 01 

-2. >26772746-01 

100 

i.ooooooooe*oo 

6,71 1037 19E-QI 

—2.412 377 96E-01 

no 

1 . lOOOOOOOE+OO 

6.47019816E— 01 

-2.40459059E-01 

120 

1.20000000F+00 

6.22989064E-01 

-2.402101 776-01 

130 

1.30000000E+00 

5.9396275 IE— 01 

-2. 4 040 2 386 E- 01 

140 

1.40000000E+00 

5. 74 899682E- 01 

-2.409752166-01 

150 

1 • 5JOOOOOOE+OC 

5. 5076376 IE-01 

-2.418680636-01 

160 

1. 600 0000 OE +00 

5.265222616-01 

-2.431149386-01 

170 

l . 70000000E+00 

5.02 14452 IE- 01 

-2.446411526-01 

180 

1.6000C000E+G0 

4.77600925E-01 

-2.464612666-01 

190 

1 • 9GOOOOOOE^OO 

4. 52 362 026 E- 01 

-2.485779136-01 

200 

7 • 0000000 OE + 00 

4.27897776E-01 

-2.510012656-01 

210 

2.100 QOOOOE + OO 

4.02676769E-01 

-2.537490436-01 

220 

2 • 200000Q0E + 00 

3. 771654 74F-01 

-2.568470336-01 

230 

2. iOOOOOOOE+OQ 

3.51327379E-01 

-2.603301356-01 

240 

2.400000 COE +00 

3.25122014E-01 

-2.642440606-01 

250 

2 • 500000 00E + 00 

2.98503763E-01 

-2.686478596-01 

250 

2. 60000000E+00 

2.71420375E— 01 

-2. 736176096-01 

270 

2. 70000000F+00 

2.4381 1036E— 01 

-2. 7925178' 6-01 

280 

2 • 80000000E+00 

2. 15603800E-01 

— 2. 8 56 7915 9F— 01 

290 

2.90000000E+00 

1. 86 71206 3E— 01 

— 2.930 70750E-01 

300 

3* OOOOOOOOE + OO 

1.57029579E— 01 

-3. 0165 82 48 E- 01 

310 

3. 1 90000 0OE *00 

1. 26 4231 76 E- 01 

-3.11763669E-01 

320 

3. 200090006+00 

9.47217U0E-02 

-3.238490556-01 

330 

3. 3000000 OE^OC 

6. 16985354E-02 

- 3. 3 86044 92 E- 01 

340 

3.400O0000E+00 

2.70421U6E-02 

-3.57114905E-01 


T AUA= 3.47422878E+00 AF ( l ) * 9 . 334 7499 OE + OO 0. 
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9TAU = 

•010 1= 10.0 

NTHETA= .30 


N 

TAU 

RHO 

RHOT 

1 

1. 00003000E-02 

9. 77512 54 IE- 01 

— 1.04 244 72 8E + 00 

13 

l .OOOOOOOOE-Ol 

9*26 830068E— 0 1 

-3. 707672 78E-01 

20 

2 . OOOOOOOOE-O 1 

8 .940416 93 E— 01 

—2. 9149 812 5E- 01 

30 

3. 00 0000 OOE— 0 1 

8 • 665243 16E- 31 

-2.5933192 IE— 01 

40 

4. OOOOOOOOE-Ol 

8.41 44974 8E— 01 

-2.419435 75E-01 

50 

5 • OOOOOOOOE-O 1 

8.17774825E— 01 

—2.3 1 26848 2 E— 01 

60 

6. OOOOOOOOE-Ol 

7.949867126-01 

-2.24247963E-01 

70 

7 • OOOOQQOOE-O 1 

7.72792333E-01 

-2.194508006-01 

80 

ti. GOOOOOOOE-Ol 

7. 51 00 726 7E- 01 

-2. 161 1 4732E— 01 

90 

9.00000000F-01 

7.29507198E-01 

—2. 13796342E-01 

100 

i.OOOOOOOOF+OO 

7.08203747E-01 

-2.12220736E-01 

110 

1.10000000E+00 

6 • 3703 1280E-01 

-2. 11209229E— 01 

120 

1*2 OOOOOOOE+OO 

6.659391 79E-01 

—2. 1 0641 477 E- 01 

130 

1.30000000E+00 

6. 44 88 705 7E— 01 

-2. 10434296E-01 

140 

1 .400000 OOE +00 

6.23841657E— 01 

-2.10529190E-01 

150 

1.50000000E+00 

6. 02 774 754E-01 

-2.10884682E-01 

160 

1 • 60 0000 OOE+OO 

5. 81 661694E— Oi 

-2. 114714526-01 

170 

1.70000000E+00 

5. 6048032 7E— 01 

—2. 122691696—01 

180 

1.80000000E+00 

5.392102 14E-01 

-2.13264392E-01 

190 

i .90000000E+00 

5.17832003E-01 

—2. 144491 83E— 01 

200 

2 • OO0OGOOOE+OO 

4.96326923E— 01 

-2.15820193E-01 

213 

2.1 OOOOOOOE+OO 

4. 74676352E— 01 

-2. 173 780 95E— 01 

220 

2 . 200000 OOE *00 

4.52 861435E— 01 

— 2. 19127276E-01 

230 

2 • 300 000 OOE+OO 

4. 308627 13E- 01 

— 2. 210 75744E-01 

240 

2 .400000 OOE +00 

4 .08 65976 IE— 01 

-2.23235222E-01 

250 

2. 50 0000 OOE +00 

3. 862308066^01 

—2. 2562 14 3 BE- 01 

260 

2.60000000E+00 

3. 635523 04E— 01 

-2.28254607E-01 

270 

2. 700000 OOE+OO 

3 • 4059846 OE— 01 

—2.311 60 1 65 E- 01 

280 

2.80000000E+00 

3. 17 340668 E- 01 

— 2. 34 36 9797 E- 01 

290 

2. 900000 OOE +00 

2.93 746 82 OE- 01 

-2.37922866E-01 

300 

3. OOOOOOOOE+OO 

2.69 7804 54 E- 01 

— 2*41 86 8384 E- 01 

310 

3.1 OOOOOOOE+OO 

2. 4539968 IE— 01 

-2. 46 267757 E- 01 

320 

3 • 2 00000 OOE +00 

2.20555772E-01 

-2.51 198646E-01 

330 

3 • 3 00000 OOE +00 

1.951913 10E-0 1 

-2. 567604 97E-01 

340 

3 .400000 OOE+OO 

1.69237668E-01 

-2.63082632E-01 

3 50 

3 .50 0000 OOE +00 

1.42611512E— 01 

-2.70336557E-01 

360 

3 • 60 0000 OOE +00 

1 • 1520981 2E- 01 

— 2. 787 54 743 E- 01 

370 

3 • 700 000 OOE +00 

8 .6902473 8E- 02 

— 2. 8866 1502 E“ 01 

380 

3. 800000 OOE+OO 

5. 752106 706-02 

— 3. 00524 804E— 01 

390 

3 . 900000 OOE +00 

2 • 68 407346 E- 02 

-3.15049329E-01 


TAUA* 3.98350950E+00 AFIi)= 1 .06660980E+01 0. 
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DTAU= .010 1= 12.0 NTHETA* .30 


N 

TAU 

RHO 

RHOT 

1 

U00000000E-02 

9. 79482068E— 01 

—9* 50677410E-01 

10 

1.00000000E-01 

9. 33308355E— 01 

-3.37379753E-01 

20 

2.00000000E-01 

9.036893 70E— 01 

-2.64921333E-01 

30 

3.00000000E-01 

8.78491927E— 01 

— 2* 3545704 IE- 01 

40 

4.00000000E-01 

8.55734784E-01 

-2* 194770 10E- 01 

50 

5.00000000E-01 

8. 34266320E— 01 

-2.09618146E-01 

60 

6.00000000E-01 

8. 13619147E— 01 

-2.03087304E-01 

70 

7.00000000E-01 

7.935265 77E— 01 

-1*98577062E-01 

80 

8* OOOOOOOOE— 01 

7. 7382 12 02E— 01 

— 1 *95390653E-01 

90 

9.00000000E-01 

7. 54390578E— 01 

— 1* 93122307E-01 

100 

UOOOOOOOOE+OO 

7. 35155168E— 01 

-1.91520002E-01 

110 

1.10000000E+00 

7. 16056312E— 01 

— 1*904 194 17 E- 01 

120 

1.20000000E+00 

6.97 049185E— 01 

— !• 89709340E-01 

130 

1 • 3 000000 OE+OO 

6. 78098445E— 01 

-1.89312306E-01 

140 

1.40000000E+00 

6. 591 75420E— 01 

—1*891731 75E-01 

150 

1 • 50000000E+00 

6.402562 16E-01 

— 1*8 *252 100E-01 

160 

1 • 60000000E+00 

6.2132041 IE— 01 

-1* 8952 003 8E— 01 

170 

1 • 70000000E+00 

6. 023501 HE— 01 

— 1* 899558 10E-01 

180 

1.80000000E+00 

5. 83329264E— 01 

- 1* 905441 20E— 01 

190 

1.90000000E+00 

5.64243133E— 01 

-1.91274204E-01 

200 

2. 000000 OOE+OO 

5.45 07789 IE— 01 

— 1*921 38902E-01 

210 

2.10000000E+00 

5.25 8202 86E-01 

— 1* 93134015E-01 

220 

2. 200000 OOE+OO 

5. 064573 72 E- 01 

-1. 942578 72E-01 

230 

2.30000000E+00 

4. 86976267E— 01 

— 1* 9551 1045E-01 

240 

2* 400000 OOE +00 

4.67363936E— 01 

— 1*96 8961 91 E- 01 

250 

2 • 50000000E+00 

4.47606988E-01 

-1.98417984E-01 

260 

2* 600000 OOE+OO 

4.27691473E— 01 

— 2* 00083141E-01 

270 

2.70000000E+00 

4, 076026 75E— 01 

—2* Q1900529E-01 

200 

2 • 80000000E+00 

3. 87 3248 86E— 01 

-2.03881360E-01 

290 

2# 900000 OOE + OO 

3.66841162E— 01 

— 2* 0603 948 2 E— 01 

300 

3.00000000E+00 

3. 46 13304 IE— 01 

— 2* 08391 785 E- 01 

310 

3* 10000000E+00 

3.25180215E— 01 

-2* 1 095 87 52 E- 01 

320 

3.20000000E+00 

3. 0396014 IE— 01 

-2* 137651 92 E- 01 

330 

3 • 300000 OOE+OO 

2. 82447565E— 01 

-2* 168412 17E- 01 

340 

3 #400 000 OOE+OO 

2. 6061393 7E-01 

-2.20223552E-01 

350 

3.50000000E+00 

2. 384266 74E-01 

-2*239573026-01 

360 

3 • 60 0000 OOE+OO 

2.15848221E— 01 

-2*280983776-01 

370 

3. 700000 OOE+OO 

1.928 348 24E- 0 1 

-2.32716863E-01 

380 

3.80000000E+00 

1.69334917E— 01 

-2* 379010 14E-01 

390 

3 • 900000 OOE+OO 

1.45 286929E-01 

—2* 437682 07E-01 

400 

4.00000000E+00 

l. 20616263E— 01 

-2 • 5046732 8E- 01 

410 

4.10000000E+00 

9.52 310032E— 02 

— 2* 58202753E-01 

420 

4* 20000000E+00 

6.90156399E— 02 

— 2* 67255892E-01 

430 

4.30000000E+00 

4. 18215542E— 02 

-2.78028639E-01 

440 

4.40000000E+00 

1.34519527E-02 

—2*911 18574E-G1 


TAUA= 4.44577991E*00 AFU> = l. 1873144 3E*01 0. 

OATE 03/07/68 TIME 19.17.02. 
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03/07/68 IRC SCOPE 3.0 6600B — 65K 02/16/68 

19.16.44. SRD4924. 

19.16.44. LRC COMPUTER COMPLEX 

19. 16. 45. JOB f 1,0644,40000. A1998, 2t 

19.16.45. MARTHA ROBI NSON , RDV122 , 1148 2011 

19.16.45. R UN (S) 

19. 16.49. SETINOF. 

19.16.50.LG0. 

19. 17.02. EOF ENCOUNTERED BY NAMELIST 

19. 17.02. EXIT 
19.17.02.SPPRINT<0UTPUT,3) 

19. 17. 03. CPU 002.206238 SEC. 

19.17.03. PPU 011.701952 SEC. 

SR04924. PRINT— PP 00806 LINES 
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XVII. APPENDIX B 

SOLUTION OF GOVERNING EQUATIONS FOR PHASE 1 


The governing equation for Phase 1 can be written as (see eq . (109) 


1 + (2p + 1 )t 2 

1(1 - p)(l + 3p) 


(B-l) 


The initial conditions for this phase are that p = 1 at t = 0. 
The motion is terminated if either the hinge circle becomes zero or the 
time, t, becomes one. The equation was solved numerically. The flow 
diagram, computer program, and numerical results follow. 

In order to avoid the singular point at p = 1, the motion of the 
hinge during the first time interval (At = 0.01) was calculated from 
the bending case. This procedure is val -1 ! as at t = 0 when p = 1.0. 
The differential equation is exactly that obtained for the pure bending 


solution. 
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SR0491 1 • IRC COMPUTER COMPLEX 

JOB* 1,0644*40000. Ai998t 2, MARTHA ROB I ^SON* RDV 122. 1148 2011 


PROGRAM PL AT ESI ( INPUT .OUTPUT! 

000003 DIMENSION CCEFSI4) ,TAU<500) , RH0(500) f RHQT <500 )• RESULT ( 2) 

CC0003 REAL IB.IE.INC 

0COOO3 COMPLEX R < 31 tTEMPC 6 ) 

00000 3 NAMEL IST/ f NPUT/OTAU* IB * IE * INC* MAXN* I PFN 

0C0CC3 PRINT 1 

000007 l FORMAT mLARGE OEFORMAT IONAL ANALYSIS OF PLASTIC*/* PLATES UNDER I 

IMPULSIVE LOADING*/* ROBINSON-KRUSZEWSK I« SRD-A1998, FEBRUARY 1968, 
2 ROV- 122* ) 

000007 CALL CAYTIM (RESULT ) 

OCOOll PRINT 2* RESULT 

0CCC17 2 FCRMAT <*0CATE*A10, 5X*T IME*A10! 

C INPUT CATA 

♦ DTAU - TIME INCREMENT 

♦ IB - INITIAL IMPULSE LOAD 

♦ IE - FINAL IMPULSE LOAD 

♦ INC - INCREMENT ON IMPULSE LOAD 

♦ MAXN - MAXIMUM NUMBER OF ITERATIONS ON VELOCITY* TIME AND RADIUS 

♦ IPFN - PRINT FREQUENCY 


OC0017 

0C0022 
000031 
000033 
000034 
000035 
OCOC40 
0C0044 
QC0064 
0C0064 
OCC 074 

0C0074 
0CCC76 
000077 
0C0100 
000103 
000107 
0C0L16 
0 CO 1 1 6 


3 READ INPUT 

CALCULATE TOTAL NUMBER OF IMPULSE LOAOS AND INITIALIZE TIME 
NN*ABS< 1 IE- IB) /INC) +1.0 
TAU( ll-OTAU 
DO 14 I *1 *NN 
N* l 

IF <1. EC.il FI* IB 
IF (I.GT.l! F I*F 1*1 NC 

IF ( < I .GT .2 ) .AND* I MODI 1*2) .NE. 0)1 PRINT 41 
41 F0RMAT<*1*) 

PRINT 4* DTAU* F I 

4 FORMAT I///* CTAU-*F5.3 t 5X*l**F5.l//) 

CALCULATE INITIAL RHO VALUE 

COEFSC 11*1.0 

COEFSI 2 ) *- l .0 

COEFSI 3)*-1.0 

COEFS<4)*l.O-CTAU/F I 

CALL FALG 1C0EF S ,3 * 0*« * TE MP» IERR) 

IF < IERR.NE.0 ) PRINT 5, I ERR 

5 FORMAT (//IOX*ERROR IN FALG SUBROUTINE*/! 

DO 6 M* 1 * 3 
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0C0120 
000124 
0C0 126 
000131 
0C0131 

OC0134 

OC0154 
OCO 167 

0CC167 
OCO 171 
OCO 173 
0C0175 
000177 


0C0220 

0C024C 

000240 

0C0244 

000247 

000251 

0CC255 

000255 

0C0261 

0C0261 

0C0264 

0C0266 


IF (AIMAG(R(M) ) ,EQ • C* 0 > GO TO 0 

6 CONTINUE 
PRINT 7 

7 FORMAT !// 10X#ALL ROOTS COMPLEX#/! 

8 RHO ( ll*REAL!RtM) ) 

CALCULATE INITIAL RHOT VALUE 

PHOT CD — !U*TAUm##2#!l.*2.*RH0!lM )/<FI#( i.-RHOIl) > •< i.#3.# 
IRHOI 1) ) ) 

PRINT 9 f N*TAU(N)# RHO! N I * RHOT I N ! 

9 FORMAT!4X*N*tliX#TAL*tlTX*RHO*»17X*RHOT#//2Xf3,3<5XEl5.8)) 

C BEGIN ITERATION ON TAUt R HU, °HOT 

10 N*N#1 

DEL RHO* RH TIN- 1) *DTAU 
RHO ! N) *RHO ( N-l ) +OELRHO 
TAU<NI*TAUCN-n#OTAU 

RHOT !NI*-!1.*TAU!N)##2#< l.#2.#RH0!N) ) l/(FI*( i.-RHO(N) )#! 1. #3.* 
lRHO(N) I ) 

C OUTPUT 

♦ RHO - RADIUS OF HINGE CIRCLE 

♦ TAU - TIME OF HINGE CIRCLE 

♦ RHOT - VELOCITY OF HINGE CIRCLE 

♦ N - ITERATION NUMBER 

IF !MOD!N« IPFN!«EQ«0) PRINT 11. N.TAU! Nl .RHOI N) , RH0TIN1 

11 FORM AT !2X 13 . 3! 5XE15 «8) I 
IF IN.GT.NMAXI GO TO 12 

IF I II .-RHGINM .LT.O.O) GO TO 14 
IF 1N.LE.2) GO TO 1C 
IF ( (TAU IN-2 )-l .0 ) *GT.0*0 ! GO TO 14 
GO TO 10 

12 PRINT 13 

13 FORMAT! //10X*MAXIMUM N REACHED#//) 

14 CONTINUE 
STOP 
END 
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PROGRAM LENGTH INCLUDING I/O BUFFERS 
OC 751 6 

FUNCTION ASSIGNMENTS 


SI/ TE ME NT ASSIGNMENTS 


1 

- CC032 3 


- 000343 

3 

- 000017 

4 

- 000353 

5 

* 000363 

7 

- 000370 

0 

- COO 13 1 

9 

- 00037A 

1C 

- 000167 

11 

- C004C5 

12 

• 000255 

13 

- 000 4*0 

14 

- 000261 

41 

* 000351 






BLOCK NAMES AND LENGTHS 


VARIABLE ASSIGNMENTS 


COEFS - 

000444 

DELRHO - 

0034 SO 

CTAU 

- 003437 

FI 

- 003445 

I 

003443 

IB 

0034C6 

It 

- 003407 

I ERR 

- 003446 

INC 

003410 

IPFN 

003441 

M 

- 003447 

NAXN 

- 003440 

N 

CC3444 

NMAX 

003451 

NN 

- 003442 

R 

- 003411 

RESULT - 

003404 

P riO 

001434 

RHOT 

- 002420 

TAU 

- U0045C 


TEMP - 003417 

START CF CCNSTANTS 
OCC 314 

START OF TEMPORARIES 
0C0415 

START OF INCIRECTS 
0C0434 

UNUSED COMPILER SPACE 
0C34O0 
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cm 


20.*t .*1, 

NQP"il 

CUNTKOL 



lk» LUi|>l* 

03**5* fki 

TiGLtS 01*202 


u.r« 

-PiJl&PiN — 

--innutis- 


--LifttLfC- 

--CONNON- 

PLiUSI 

COOICO 




MSTFW 

0076 I ft 




cctpft 

OlCSfel 




( **T |M 

010716 




IMCTN 

L 1 C 7*7 




1 U f. 

Pl^ll 




f 1 CN 

CH?7| 




(1 TP* 

01*11* 




m ni t* 

01*1 s* 




CPC 

C 1 5 J? 7 




NT»y — 

— iC 0 «f ss- 




PI AttSl 

COC 101 




i*KT»V 

CCI617 

PLiTCSl 

00010? 


MMF* 

CIOCC? 

GLlPTC 

010575 




INPUTN 

011*20 

011*30 



KGOEi 

0151*0 


ST!f#«t 

CO 7 7*7 




jTSftPP 

CO 7 7 7S 




ISC 

001671 

PLiTtSl 

000367 




FilG 

013120 


M(> 

CC772? 

PLiTtSl 

000)65 


MIT 

GC 7 71 * 

INPUTN 

012207 



ot.ru/ 

OLTPTC 

01C676 




1 ftPUTf 

011*31 




RJOtP 

0151*1 


CO IP tl 

01C363 

PLiTtSl 

000105 

000106 




000166 

000170 




000227 

000230 




000266 

000)26 




000)67 

000)60 

tiVtl" 

010717 

PLiTtSl 

000110 


ISPuTN 

010760 

PLilESl 

000121 


MiG 

01221* 

PLiTtSl 

C00206 


ffft*PMi, 

oim* 





Cl 3675 




PQ6PIL. 

CIU1) 

CL T PTC 

01061) 


•DPPU. 

C 1 36* * 




lit. 

0136*5 

OLTPTC 

010616 

01C6I5 

ctot. 

C 1 16 31 




(PIN. 

01)275 

S VST E* 

C10257 




OLTPTC 

01060) 




INPUTN 

011)12 


!IU. 

Cl 3*07 

CUTPTC 

01Q650 


Ct 1*6 

01*115 

OUT PTC 

010571 




INPUTN 

010756 


*Ui>t * 

C1*M« 

CUTPTC 

010565 

010*21 

CK 

016 17 6 

DiVTH 

01072* 

0107)0 


MC02 US«.*1 

(PCM Cl 4 32 7 

(Hf* CIS 1 * ft 

IPC 1 *** C|*I^ 

UNSiMSHtU MT fcMNAL S- 


0001 PO 015566 000030 000009 

♦ CUl f M loan- -III A l 060 — Si N6 C.CNN-- (IMIh.. 


PIFIMNCF* 


00011) 

000116 

000116 

0001*2 

00 01*7 

OOC172 

COCITt 

QJ0212 

00021* 

0002 l 5 

000255 

000257 

000261 

00026) 

0002 6 5 

00030 

0001)2 

0001)* 

000)36 

000117 


01 C6*i 


OIOS* 


•if FI *t NC I S 
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LARGE DEFCRMATICNAL ANALYSIS OF 
PLATES UNCEP IMPULSIVE LOADING 
RCBlNSCN-KPLSZEkSKl . SRD-A1998, 


PLAST IC 


FEERUARY 1968, RDV-I22 


DATE C3/G6/68 


TIME 20.41.43. 


D TAll= .CIC 


t= 2.0 


N 

TAU 

RHQ 

RHOT 

1 

1.00CC0CC0E-02 

9. 4935465 OE- 01 

-2 *5663 3899E+00 

10 

l.OOOCOOCGE-Ol 

6. 2 6 SC3 2426-0 1 

-8 • 60407 945E -0 1 

20 

2*0CCC00C0E-01 

7 * 52748 29 3 E-0 1 

-6 *82851812 E- 01 

30 

3 .OOOCOCCOE-Ol 

6.67216557E-01 

— 6* 33697 1 90E— 0 1 

40 

4.000090COE-01 

6.24269401E— 01 

-6*29870 165E-0 1 

5C 

5.QOCCOCCOE-G1 

5.6C554 158E-01 

—6 *492 7 78 74 E-0 1 

60 

6 *0 00 0000 OE -01 

4.94193799E-01 

-6*832 13 66 9E -01 

7C 

7.CCCC0CCCE-01 

4*2 39445 39E— 0 1 

— 7*2799 8205 E- 01 

80 

8. )QGCOOCCt-Ci 

3.48 763 53 3E— C 1 

— 7* 82825 86 C E-0 1 

9G 

S.COOOOOOOE-Oi 

2 *67 592 1 C4E— 0 1 

— 8*49573 7C9E-01 

100 

l.CGCCOCCCE+OO 

U79C1262IE-01 

—9 *34323791 E— 0 1 

TA 0= 

.01C 1= 4.0 



N 

TAU 

RHO 

RHOT 

L 

l • CGCCQGC CE-02 

9 *64 325 054 E—0 I 

-1.30062C63E+00 

1C 

l.COOOOOCOE-Ol 

8.8C6C7841E-01 

- 5. 90846328E -0 1 

20 

2 .OOOOOOC OE-0 1 

3.28619291E-01 

-4.62953323E-01 

30 

3.000CGCC0E-C1 

7* 943798 86 E— 01 

— 4*2 571 955 8 E-0 1 

40 

4 .00000000 E-0 1 

7 *42234 1 59E— 01 

-4 . 2006082 3E-Q 1 

5C 

5.CCCCCCCCE-01 

6.9« -7598oE-01 

-4.299647436-0 1 

60 

6 .OOOOOOOuE-Ol 

6 • 5t>C90 24 3E— 0 1 

— 4. 48 7544C4E -0 1 

7C 

7.CCCC0CC0E-01 

6.1015748CE-G1 

-4 • 7 3056 340 E-0 1 

80 

b *000000006-0 1 

5*6 1618031 E— 01 

— 5.010 38 2 04 E-0 1 

90 

S*Q0GC0CCCt-01 

5 • 10152 76 l E-0 1 

-5*317 39944 E-0 1 

100 

1 *000 0000 OE +00 

4.555 10761E-01 

-5*647856986-01 
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QT AU= .01C 1= 6.0 


N 

TAU 

RHO 

RHOT 

1 

1 .000000006-02 

9 . 7C 520 305E-0 1 

-1.46522169E+00 

10 

i.COCCOOOOE-Gl 

5.03C65367E-0 I 

-4.765505756-01 

20 

2.00000000E-01 

8.61 227 120E-0 l 

-3* 7 1626 28 8E- 01 

30 

3.C00CCCC0E-01 

8 . 25 765 008 E-0 1 

-3 *40734 5136-0 1 

AO 

4.000C00006-01 

7.92C63572E— 01 

-3*355635866-01 

50 

5 •OOOOOOOOE-O 1 

7 .58 248952 E-0 1 

- 3 .4297062 1E-0 1 

6C 

fc.GOOOOOCOE-Cl 

7.2 3345 1 59E-01 

-3.574299456-01 

70 

7 •OOGOOOOOE-O 1 

6 .66786 77 IE-01 

-3.761015936-01 

80 

8.00000CC06-01 

6.482379G1E-01 

-3.97382427E-0 1 

90 

9*000000005-01 

6.C7477632E-01 

-4. 203439486-01 

IOC 

1 .CCGCQCOOE+OO 

5.64 3658 18E-0 1 

-4 .444706C 5 E— 0 1 


DTAU= .01C 1= 8.0 


N 

TAU 

RHO 

RHOT 

1 

1.00000000E-02 

5.74E4125SE-01 

-1. 2663746 9E •♦•00 

10 

1.C0C00000E-01 

9. 16329 031 E- 01 

-4.09781482E-01 

20 

2.0C000CC0E-01 

8. 8C 397 088E-0 1 

-3.18725378E— 01 

30 

3.00000000E-01 

8.50C05267E-01 

~2.9 17932C7E— 0 i 

40 

4.CCCCCCC0E-01 

8.21156 722E-01 

-2 . 8 7 1 1 8 1 76 E-0 1 

50 

5.000000C0E-01 

7.92 231636E— 0 1 

-2. 9329128 5E -01 

60 

6.O0CCOO0OE-O1 

7 .62 39044 8 E-0 1 

-3.05500383E-01 

70 

7.COCOOCCOE-01 

7.21 1536C6E-01 

-3.21258145E-01 

8C 

e.cocooocoE-oi 

6. 982 36 316 E— 0 1 

-3.39145222E-01 

90 

5.COCCOOCOE-01 

6. 63467C56E-01 

-3. 5832015 IE— 01 

100 

1.00000000E+00 

6. 26741673 E-01 

-3. 7826809 8E-01 
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0 T A U= . C 1 C 1= 10.0 


N 

TAU 

RHO 

RHCT 

l 

1.C00C0OC06-02 

9*7 75 1254 IE— 01 

-1.131136596+00 

10 

l .CQ0000C0E-01 

9*252290656-01 

- 3 *6477492 IE- 0 1 

20 

2.0CCC0CC0E-01 

6* 93 3692 78 E— 01 

-2*832400566—0 1 

30 

3.00000000E-01 

8*66373491 E— 01 

-2.590661326-01 

40 

4 *0000000 0E -01 

8*40766 777 6-01 

-2. 547920946-01 

50 

5.000C0C0CE-01 

8. 15 101 7 1 5E- 01 

-2.601993976-01 

60 

6*000000C0E-0l 

7*886301396-01 

-2*709 705 10E-0 1 

70 

7.CCCCG0C06-01 

7*608272106—01 

—2*84867606 E-0 1 

80 

8 *00000 OOOE-0 1 

7.31 743604E-0 1 

-2.0C6C7661E-01 

90 

8.CCCC0QG0E-Q1 

7.00932 7616-01 

-3*174229156-01 

100 

l.OOOOOCCOE+OO 

6. 664098 74E-01 

-3.348524626-01 


0TAU= .C1C 1= 12.0 


N 

TAU 

RHO 

RHOT 

1 

i.ooccoocoe-G 2 

9. 7S482 068 E-0 1 

-1 *031546326+00 

10 

1 .OOCOOOOOE— 01 

9.3 1946283E— 01 

-3.318349246-01 

20 

2. CCO COO COE— 01 

9. C 2 889 25 7 E— 01 

-2.573526226-01 

30 

3 .OOOOOOOOE— 0 1 

8. 7E3686C1E— 01 

-2 * 3523 8060E-0 1 

40 

4.0CCC0000E-01 

8.55 12078 IE— 01 

-2*312356496-01 

50 

5. 000 COCO OE— 01 

6.3 1E25262E— 01 

-2. 361586836-01 

6C 

6.00000000E-01 

8 ,C 7800665 E— 01 

-2.45907867E-01 

70 

7.0CCC0CC0E-01 

7.82661535E— 01 

-2.58482927E-01 

80 

8 .OOOOOOOOE— 01 

7.56 183372E— 01 

-2.72705358E-01 

9C 

9.00000CCOE-01 

7.2 8236 22 7 E -01 

-2.878666096-01 

100 

l.OOCCOOCOE+OO 

6.9E747338E— 01 

-3.035370766-01 
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DTAU= .CIC 1= 14.0 


N 

TAli 

RHO 

RHOT 

1 

l .OOGCOOCOE— 02 

9.81011420E-01 

— 9.54283778E-01 

10 

1 .OOCCOOOOE— 01 

9.37C73938E-01 

—3. 06396242 E-01 

20 

2.00000000E-01 

9. 10256 151 E— 01 

—2. 374073 13E-0 1 

30 

3.O0OC00CCE-01 

8.87641213E— 01 

-2.16904737E-01 

40 

4.00000000E-01 

8 .66 2C7 554E— 0 1 

-2. 13214769E-01 

50 

5 .COCCOOOOE— 01 

8. 44733200 E-01 

-2.17687642E-0 1 

60 

6.00000000E-C1 

8. 22588230E— 01 

-2.26661812E-01 

70 

7.0C000C00E-01 

7. 99417 35 2E— 01 

-2.38234367E-0 1 

80 

e.ooooocccE-ci 

1.75C14 706 E-01 

-2.513C9793E-01 

90 

9 .00000000E— 0 1 

7 .492624 1 3E-0 1 

-2.65227619E-01 

IOC 

i.cocooccoe*oo 

7.22C960C1E-01 

—2.79585 154 E-0 1 


0T AU= .QIC 1= 16.0 


N 

TAU 

RHO 

RHOT 

L 

l.OOOOOOOOE— 02 

9.82 243 330 E— Cl 

— 8.92092324E-01 

10 

1 . COOCOOCOE— C 1 

9.41197612 E— 0 1 

—2 .85992394E-0 1 

20 

2.00000000E— 01 

9. 16 174345E— 01 

— 2. 2 143825 3E -01 

30 

3.CCOOOOCOE— 01 

8.95084436E-01 

—2.02241 39 5 E-01 

4G 

4.0COCOOCOE— Cl 

8. 7 5 10143 IE— 01 

— 1 • 98769 786 E-0 1 

50 

5.00000000E— 01 

8. 55082496 E-01 

—2.0 292833 7E- 01 

60 

6 • OCOCCCOOE— 01 

8. 34439220E— 01 

—2. 11288325 E-01 

70 

7.00000000E— 01 

8.12640285E— 01 

—2. 2206645 7E-01 

80 

e.ocoooocoE-oi 

7.9CC94500E-01 

—2.34235 15 2 E-0 1 

90 

9.CCOCOOOOE— 01 

7.66093258E-C1 

-2.47173616E-01 

100 

i.coccoocce*oo 

7. 4C 778 205 E-0 1 

—2.60502 145 E-01 
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DTAU= .C1C 1= 18.0 


N 

TAU 

RHO 

RHOT 

1 

1 .COCCOOCOE— 02 

9.82 263156E— Cl 

— 8.40637554E-01 

10 

1 .00000000E-01 

9 .4 4 6066 69E— 01 

-2.69159456E-01 

20 

2 .000000 COE- Cl 

9.21C63 053E-01 

-2.08282849E-01 

30 

3.000000C0E-01 

9.01228986E— 01 

— 1. 9017121 0E-01 

40 

4.00G00000E-01 

8. 82439802 E-01 

-1.86884508E-01 

50 

5.00000000E-01 

e.63fcl8262E— 01 

— 1.90 787 59 7E— 01 

60 

6.CC0C0000E— 01 

8 .44 2 10 l 32 E— 01 

— 1.98645 14 IE- 01 

70 

7 .OGCC'OOOOE— Cl 

8.2290J619E— Cl 

—2 • 08773523E— 0 l 

80 

8 .OOOOOOOOE-O 1 

8.02520068E-01 

-2,20201 743E-01 

90 

9 .COCCOOCOE— 0 l 

7. 79 9 576 82 E— 01 

—2.3 2342675 E-0 1 

100 

1.000000C0E*00 

7.56 163042E— 01 

— 2. 448 3626 7E— 01 


DTAL= .C1C 1= 20.0 


N 

TAU 

RHO 

RHOT 

1 

l.OCCCCOCOE -02 

9. E4125486E— 01 

—7 .97150195E— 0 1 

10 

1 .OOOOOOC OE-O l 

9.41486469E-01 

-2.5496681 3E-01 

20 

2-COOOOOCOE-Ol 

9.25189353E-01 

— 1.9720 3850 E- 01 

30 

3.OOOC0OCOE-01 

9. 06 4 12 5 7 3E- 01 

— 1. 8001268 OE— 01 

40 

4. OOOOOOOOE-O 1 

8.6E62T966E— 01 

— 1.7 68 8495 OE- Cl 

50 

5.CCCC0CC0E— 01 

8. 7C 8 13 749 E— 01 

— L.80575023E-0 1 

60 

6.00000000E— 01 

8.52444521E-01 

—1. 8801 151 9E— 01 

70 

7.0CCOOOOOE-01 

8.332253C1E-01 

-1.97595440E— 01 

80 

8.00CC00C0E— 01 

€. 12986759E— 01 

—2 • 08404C89E-0 1 

9C 

9. OOOOOOOOE-O 1 

7 .91633808E— 0 1 

— 2.198 79C69E— 01 

100 

l.OCCCCCCOE+OO 

7.69116623E-01 

—2.31677 393E— 0 1 



02/C6/68 LRC SCOPE 3.0 6600B — 65* C2/16/68 

2C.41 .33. SRC491 1 . 

20.41.33. LRC COMPUTER COMPLEX 

20.41.33. JOB, 1,0644.40000. A1998. 2, 

20.41.33. MARTHA RCBINSGN.RDV122. 1148 2011 

20.41 .34. RUN( S ) 

2C.41.38. SETINOF. 

20.41.35. LGC. 

20. 41. 45. STOP 

2 C. 4 1.4 5. SPPRl NT! OUTPUT. 3» 

20. 41. 47. CPU 000.950224 SEC. 

20. 41. 47. PPL C06.35640C SEC. 

SRD4911. PRINT-PP GO 392 LINES 
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XVIII. APPENDIX C 
VALUES OF F(n) 

The value of the function F(n) is needed to evaluate the final 
deformation. Since the argument of F(n) is limited between 0 and 1, 
it was expedient to generate a table for this function. 
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X 

F (X) 



0.00 

0 . 



.01 

4.02093263E-02 

.51 

4. 96479867E+00 

.02 

8. 08746017E-02 

.52 

5.24893643E+00 

.03 

1.22051904E-01 

.53 

5. 55417411E+00 

.04 

1.63797750E-01 

.54 

5. 882 70543E+00 

.05 

2.06169417E-01 

.55 

6 . 2370102 7E+00 

.06 

2. 49225241E-01 

.56 

6.61989978E+00 

.07 

2 . 930 2492 5E-0 1 

.57 

7.03457006E+00 

.08 

3.37629834E-01 

.58 

7.48466593E+00 

.09 

3.83103291E-01 

.59 

7.97435731E+00 

.10 

4.29510892E-01 

.60 

8.50843084E+00 

.11 

4.76920811E-01 

.61 

9.09240044E+00 

.12 

5.25404131E-01 

.62 

9. 73 264126E+00 

.13 

5.75035184E-01 

.63 

1. 04365528E+01 

.14 

6.25891907E-01 

.64 

1.12127588E+01 

.15 

6.78056228E-01 

.65 

1.20713527E+01 

.16 

7.31614470E-01 

.66 

1 . 302 42020E+01 

.17 

7.86657786E-01 

.67 

1.40853277E+01 

.18 

8.43282629E-01 

.68 

1.52713783E+01 

.19 

9.0159 1261E-01 

.69 

1. 660223 01E+01 

.20 

9.616923 00E-01 

.70 

1.81017486E+01 

.21 

1.02370132E+00 

.71 

1.97987645E+01 

.22 

1.08774151E+00 

.72 

2. 17283338E+01 

.23 

1.15394437E+00 

.73 

2. 39333795E+01 

.24 

1 . 22245053E+00 

.74 

2.64668506E+01 

.25 

1.29341057E+00 

.75 

2.93945880E+01 

.26 

1, 36698601E+00 

.76 

3.27991738E+01 

.27 

1.44335031E+00 

.77 

3.67851575E+01 

.28 

1.52269007E+00 

.78 

4.14862437E+01 

.29 

1.60520631E+00 

.79 

4. 70753123E+01 

.30 

1. 69111587E+00 

.80 

5.37785952E+01 

.31 

1.78065300E+00 

.81 

6.18960611E+01 

.32 

1.87407114E+00 

.82 

7.18312529E+01 

.33 

1.97164490E+00 

.83 

8.41358234E+01 

.34 

2.07367220E+00 

.84 

9.95774710E+01 

.35 

2 . 18047681E+ 00 

.85 

1.19246106E+02 

.36 

2.29241109E+00 

.86 

1.44724337E+02 

.37 

2. 4098590 8E + 00 

.87 

1.78369766E+02 

.38 

2.53324005E+00 

.88 

2.23799110E+02 

.39 

2. 66301240E+00 

.89 

2 . 86752518E+02 

.40 

2. 79967822E+00 

.90 

3 . 76710856E+02 

.41 

2. 9437882 8E + 00 

.91 

5.10094046E+02 

.42 

3.09594788E+00 

.92 

7.17020114E+02 

.43 

3 , 256823 32E+00 

.93 

1.05679676E+03 

.44 

3. 4271495 1E+00 

.94 

1.657233 04E+03 

.45 

3.60773845E+00 

.95 

2. 82847838E+03 

.46 

3.79948909E+00 

.96 

5.45742996E+03 

.47 

4.00339864E+00 

.97 

1.27818794E+04 

.48 

4. 22057561E+00 

.98 

4.26332497E+04 

.49 

4.45225484E+00 

.99 

3.37140837E+05 

.50 

4. 69981495E+00 





